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Abstract: A comprehensive treatment of black hole thermodynamics in two-dimensional 
£^ dilaton gravity is presented. We derive an improved action for these theories and construct 

the Euclidean path integral. An essentially unique boundary counterterm renders the im- 
proved action finite on-shell, and its variational properties guarantee that the path integral 
has a well-defined semi-classical limit. We give a detailed discussion of the canonical en- 
semble described by the Euclidean partition function, and examine various issues related to 
stability. Numerous examples are provided, including black hole backgrounds that appear 
in two dimensional solutions of string theory. We show that the Exact String Black Hole 
is one of the rare cases that admits a consistent thermodynamics without the need for an 
external thermal reservoir. Our approach can also be applied to certain higher-dimensional 
5— i black holes, such as Schwarzschild-AdS, Reissner-Nordstrom, and BTZ. 
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1. Introduction 

Two-dimensional gravity provides a useful expedient for testing ideas about quantum grav- 
ity in higher dimensions. Technical simplifications in two dimensions often lead to exact 
results, and it is hoped that this might help to address some of the conceptual problems 
posed by quantum gravity. At the same time, these simplifications remove many of the fea- 
tures that make gravity interesting. Striking a balance between models that are tractable 
and models that seem relevant is an art in its own right. 

In this paper we study black hole (BH) thermodynamics in two-dimensional models 
of dilaton gravity. In the absence of matter these models can be solved exactly, but they 
contain no dynamical degrees of freedom. This seems too restrictive, so we sacrifice inte- 
grability and consider models that are coupled to matter. A detailed description of the 
matter is not necessary for our purposes. We only require the existence of some propagat- 
ing degrees of freedom, analogous to the radiation emitted by a BH in higher dimensions. 
We then perform a semi-classical analysis and obtain a comprehensive description of BH 
thermodynamics for arbitrary models. Despite the fact that we work in two dimensions, 
these systems exhibit a wide range of interesting thermodynamic phenomena. To some 
extent this is due to our consideration of 'quasilocal' thermodynamic potentials describing 
BHs in a finite cavity. Such a treatment reveals aspects of BH thermodynamics that are 
not visible in an asymptotic analysis. Several universal properties of dilaton gravity BHs 
are identified, as well as structures that are common to families of models. Our approach 
allows us to study the thermodynamics of backgrounds, both approximate [1] and exact [2], 
associated with two-dimensional solutions of string theory. We are also able to extend our 
analysis to spherically symmetric BHs in higher dimensions. In particular, we obtain an 
expression that describes the free energy of both asymptotically flat and asymptotically 
Anti-de Sitter (AdS) BHs. 

Dilaton gravity in two dimensions is conventionally described by the (Euclidean) action 

/=-— L- [d 2 x^ \XR-U(X) (VX) 2 -2V(X) ] -— L- fdx^jXK. (1.1) 

We define the dilaton X in terms of its coupling to the Ricci scalar by the expression X R. 
Different models are distinguished by the kinetic and potential functions U(X) and V(X), 
cf. e.g. [3] for a review. The remaining part of the action is the analog of the Gibbons- 
Hawking- York boundary term [4, 5] , where 7 a b is the induced metric on the boundary and 
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K is the trace of the extrinsic curvature. Our conventions for curvature tensors and other 
geometric quantities can be found in appendix A. To study the thermodynamics of these 
theories we employ the Euclidean path integral, 

Z = J 9g®Xexp(~I\g,X] \ . (1.2) 

The path integral is evaluated by imposing boundary conditions on the fields and then 
performing the weighted sum over all relevant space-times ( A4 , g) and dilaton configurations 
X. In the semi-classical limit it is dominated by contributions from stationary points of 
the action. This can be verified by expanding it around a classical solution 

I[g d + Sg, X d + 8X) = I[g ch X cl ] + SI[g d ,X d ; Sg, SX] + ~Pl\g d ,X d ; Sg, 8X) + ... (1.3) 

where 51 and 5 2 I are the linear and quadratic terms in the Taylor expansion. The linear 
term is assumed to vanish on-shell. Then, if the leading term is finite and the quadratic 
term is positive definite, the saddle point approximation of the path integral is given by 

Z ~ exp (-~I[g d ,X d ] \ J MgMX exp (-~8 2 I[g cl ,X d ;8g,5X]\ . (1.4) 

This is interpreted as the partition function in the canonical ensemble, with the temperature 
set by the path integral boundary conditions. 

The transition from (1.2) to (1.4) is complicated by the fact that the action (1.1) 
does not possess any of the properties listed above. First, the on-shell action diverges. 
This is familiar from calculations in higher dimensions, and is commonly addressed using 
a technique known as 'background subtraction' [5,6]. Second, the linear term in (1.3) 
does not vanish for all field configurations that contribute to the path integral [7, 8] . This 
is due to boundary terms whose contribution to 51 are usually not considered in detail. 
Schematically, these boundary terms take the form 
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on— shell 



I dx ^7 U ab hah + vrx 5X 1 + . (1.5) 

JdM 1 J 



Even when we consider models where the boundary conditions imply 5j a b — ► and 5X — ► 
at dA4, the coefficients of these terms tend to diverge so rapidly that 51 does not vanish. 
Finally, the Gaussian integral in (1.4) is often divergent. In that case the canonical partition 
function is not well-defined. Instead of describing the thermodynamics of a stable system, 
(1.4) yields information about decay rates between various field configurations with the 
specified boundary conditions [9]. 

We solve the first two problems by replacing (1.1) with an 'improved' action T. The 
improved action is related to (1.1) by a boundary counterterm. The counterterm, which 
is essentially unique, is obtained by extending the results of [10] to arbitrary models of 
dilaton gravity. The improved action is finite on-shell, and 5T vanishes for any variation 
of the fields that preserve the path integral boundary conditions. However, it does not 
address the remaining problem. The quadratic term 5 2 T may result in a Gaussian integral, 
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as in (1.4), that is divergent. This is directly related to the thermodynamic stability of the 
system. Essentially, the density of states grows so rapidly that the canonical ensemble is 
not defined. We deal with this in the same manner as York [11]. The BH is placed inside 
a cavity that couples the system to a thermal reservoir, with boundary conditions fixed at 
the wall of the cavity. Following Gibbons and Perry [12] , we identify the location of the 
wall with the value of a conserved dilaton charge. The canonical ensemble that emerges 
from this procedure is well defined if and only if the specific heat of the system is positive. 
Typically this requires a finite cavity, but there are some notable exceptions. Having 
resolved the three problems described above, the approximation (1.4), with T replacing /, 
gives the canonical partition function for the system. 

We proceed as follows: in section 2 we review BH solutions of two-dimensional dilaton 
gravity [13, 14], and examine the action (1.1) and its variation evaluated on these back- 
grounds. We then derive the improved action T and verify its properties. The Euclidean 
path integral is used to define the canonical partition function in section 3, from which all 
of our results on thermodynamics follow. We obtain expressions for thermodynamic quan- 
tities in an arbitrary model of two-dimensional dilaton gravity, establish universal results 
like the first law of BH thermodynamics, and examine several issues related to stability. In 
section 4 these results are applied to two significant families of models. Section 5 addresses 
BHs that appear in models obtained from string theory. This includes a treatment of the 
Exact String BH [2]. The results obtained in section 3 can also be applied, in certain 
cases, to BH solutions of higher dimensional theories. This is discussed in section 6, where 
we study spherically symmetric BHs with or without a cosmological constant, and the 
toroidal Kaluza-Klein reduction of the BTZ BH [15] . Finally, in section 7 we summarize 
our results, discuss further applications, and suggest possible generalizations. A collection 
of appendices contain our conventions, a review of various models, and additional technical 
results. 

Before proceeding to section 2, we point out a few important conventions. Euclidean 
signature is used throughout this paper. Nevertheless, we frequently use terms such as 
'Minkowski' and 'horizon' as if we were working in Lorentzian signature. The dimensionless 
Newton's constant that appears in the action is set to 8ttG2 = 1. It is restored in certain 
expressions, when necessary. 

2. Semi-Classical Approximation of the Dilaton Gravity Path Integral 

In this section we consider the Euclidean path integral formulation of dilaton gravity in 
two dimensions. Our analysis shows that the standard form of the action is not stationary 
under the full class of field variations that preserve the boundary conditions. We construct 
an improved action that solves this problem, and obtain the semi-classical approximation 
of the corresponding path integral. 

2.1 Equations of Motion 

The equations of motion are obtained by extremizing the action (1.1) with respect to small 
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variations of the fields that vanish at the boundary. This yields 

U{X) V^XV U X - l - V [/(1)(VJ) 2 - g^V{X) + V M V,X - g^V 2 X = (2.1) 
R + d x U{X)(VX) 2 + 2 U{X)V 2 X - 2 d x V{X) = . (2.2) 

Solutions of these equations always possess at least one Killing vector whose orbits are 
curves of constant X [16, 17]. We fix the gauge so that the metric is diagonal. The 
solutions I 

X = X(r) , ds 2 = £(r) dr 2 + — dr 2 , (2.3) 

with 

d r X = e- Q{x) (2.4) 

are expressed in terms of two model-dependent functions, 

Q{X) := Q + fdXU{X), (2.6) 



w(X) :=w -2 dX V(X) e Q(x) . (2.7) 

Here Qq and wq are constants, and the integrals are evaluated at X. Notice that wq 
and the integration constant M contribute to £(X) in the same manner. Together, they 
represent a single parameter that has been partially incorporated into the definition of 
w(X). By definition they transform as wq — * e^°wo and M — * e A ®°M under the shift 
Qo —> Qo + ^Qo- This ensures that the functions (2.4) and (2.5) transform homogeneously, 
allowing Qo to be absorbed into a rescaling of the coordinates. Therefore, the solution is 
parameterized by a single constant of integration. With an appropriate choice of wo we 
can restrict M to take values in the range M > for physical solutions. As evident from 
(2.3) the Killing vector described above is d T . For brevity we r efer to y/£(X) from (2.5) 
as the "Killing norm" . If it vanishes we encounter a Killing horizon. The square of the 
Killing norm for the solution with M = is 

Co(X) := w{X) e Q W . (2.8) 

This solution, which we often refer to as the 'ground state', plays an important role through- 
out this paper . 

2.2 Black Holes 

The constant of integration that labels solutions of the equations of motion is analogous 
to the mass parameter that appears in the Schwarzschild metric. This constant has been 



1 If the function V(X) has one or more zeroes at finite values of X there is a second, inequivalent family 
of solutions with constant dilaton. In that case, the identification of the actual ground state of the model 
becomes more involved. This issue will be addressed in subsection 3.11. 
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partially absorbed into the definition of w(X) in such a way that M = represents the 
ground state of a particular model. We refer to solutions with M > 0, which may exhibit 
horizons, as BHs. 

In the models we consider the metric function is no n- negative over a semi-infinite 
interval 

X h < X < oo . (2.9) 

The upper end of this interval corresponds to the asymptotic region of the space-time. The 
function w{X) generally diverges in this limit 

lim w(X) -» oo (2.10) 

so that the asymptotic behavior of the metric is typified by the ground state solution (2.8). 
The lower end of the interval (2.9) is either the origin, or a Killing horizon. Assuming that 
the function e^ x ^ is non-zero for finite values of X, the location of the horizon is related 
to the parameter M by 

w{X h ) = 2M . (2.11) 

If this condition admits multiple solutions then is taken to be the outermost horizon, 
so that w(X) > 2M for X > X h . 

Regularity of the metric at the horizon fixes the periodicity r ~ r + (3 of the Euclidean 
time, which is given by 

4tt 
w'{X) 
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(2.12) 

x h 



I'h 

The inverse periodicity is related to the surface gravity of the BH by 2-nj3~ 1 = k. In an 
asymptotically flat space-time, j3~ l is also the temperature measured by an observer at 
infinity. In a slight abuse of notation, we denote this quantity by T 
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(2.13) 



This should not be confused with the proper local temperature, which is related to (3 1 by 
a position dependent red-shift (the 'Tolman factor' [18]) 

T{x)= vm rl - {2M) 

The proper temperature at infinity coincides with (2.13) only if £(X) — ► 1 as X — > oo. 

At this point it is instructive to evaluate the action for the BH solution, which we 
expect is related to the semi-classical approximation of the BH partition function. The 
bulk terms in (1.1) are straightforward, but we need to specify how the boundary is treated 
in order to calculate the Gibbons-Hawking- York term. It is apparent from (2.3) that the 
only component of dM that contributes to this calculation is a surface orthogonal to the 
coordinate r. Depending on the choice of coordinate system this may occur either at some 
finite value of r, or in a limit such as r — > oo. In either case, the boundary is assumed 
to correspond to X — > oo, which means that the boundary terms in the action should be 
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calculated via a suitable limiting procedure. We implement this by placing a regulator 
X < X Teg on the dilaton, and treating the isosurface X = X reg as the boundary. Then the 
on-shell action is 

I reg =p(2M- w{X veg ) - 2 irX h r) . (2.15) 

The regulator should be removed by taking the limit X TCg — ► oo, which recovers the full 
BH space-time. At this point we encounter a problem: the asymptotic behavior (2.10) of 
w(X mg ) causes the action to diverge. How should this term be addressed? One possibility 
is that it should be discarded, and the finite remainder then interpreted as the action for the 
space-time. This strategy is unsatisfactory for two reasons. First, it fails to reproduce the 
correct thermodynamics of the BH. Second, as we shall demonstrate in the next subsection, 
the variational properties of the action (1.1) are not consistent with our assumptions about 
the path integral. Neither issue is resolved by the sort of ad-hoc subtraction described 
above. 

2.3 Variational Properties of the Action 

The equations of motion were described in section 2.1 as the conditions for extremizing 
the action. In this section we show that this is true for 'localized' variations of the fields, 
but not for generic variations that preserve the boundary conditions. This has important 
consequences for the path integral formulation of the theory. We begin by considering 
small, independent variations of the metric and dilaton. The corresponding change in the 
action is 

SI = -\ I d 2 Xy /g\ 5g^ + £ x 5x} + [ dx ^ \ ir ab 5 lab + ir x 5X ] . (2.16) 
2 Jm 1 1 JdM 1 1 

The terms £^ v and Ex m the bulk integral are the left-hand-sides of (2.1) and (2.2), 
respectively. In addition to the bulk term, there is also a boundary term that depends on 
momenta conjugate to 7 & and X 

K ab = ~\ l ab n^V^X , ttx = U(X) n^V^X - K . (2.17) 

These momenta are defined with respect to evolution of the fields along the direction normal 
to the boundary, indicated here by the vector n^. When evaluated on a solution of the 
classical equations of motion, the change in the action (2.16) should vanish for any variation 
of the fields that preserves the boundary conditions. A direct calculation, which we shall 
carry out below, shows that this is not the case. Instead, we find that 51 ^ unless we 
restrict our attention to localized variations of the fields. In this context a field variation 
is 'localized' if its asymptotic behavior satisfies some condition that is more restrictive 
than simply preserving the boundary conditions. Essentially, the change in the action only 
vanishes for variations that fall off more rapidly at infinity than required by the boundary 
conditions. 

The claim that 51 does not vanish for physically reasonable variations of the fields 
may seem odd. The Gibbons-Hawking- York term in (1.1) is usually motivated by stating 
that it leads to an action with a "well defined" variational principle. But the Gibbons- 
Hawking- York term only ensures that fields obey Dirichlet conditions at dM. It does 
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not guarantee that the boundary term in SI vanishes for arbitrary 5^f a b an d SX that 
preserve these boundary conditions. This issue has been studied in higher dimensions for 
gravitational theories with both asymptotically fiat [7, 8] and asymptotically AdS [19, 20] 
boundary conditions. It is especially important for the path integral formulation of the 
theory. Given a solution (£ c ;, X c {) of the equations of motion, the path integral also receives 
contributions from field configurations whose asymptotic behavior is (£ c / + 5£,X c i + SX), 
where S£ and SX are generic variations consistent with the boundary conditions. The 
h — > limit of the path integral is dominated by solutions of the classical equations of 
motion if, and only if, SI vanishes for all such variations. 

We show now that the action (1.1) does not satisfy this condition. Evaluating (2.16) 
for a solution of the form (2.3) gives 



SI TC g = I dr 



■ l - 8 r X 5£ + (u(X) £(X) d r X - l - d r i 



(2.18) 



The subscript on SI indicates that the boundary term must be evaluated using a regulator, 
as in the previous subsection. The boundary term should vanish as the regulator is removed 
to infinity. To show that this is not the case, it is enough to set SX = and consider the 
variation Using (2.4), the coefficient of 5£ is proportional to e~®^ x \ Depending on 
the model, this term may either diverge, approach a constant, or vanish at dM. Now we 
need to specify how a generic variation 5£ is allowed to behave at dM. The models we 
consider allow for a wide range of asymptotic behavior. In particular, £ might diverge in 
the X — > oo limit. Thus, we cannot assume that <5£ vanishes at dM, as we would if £ 
approached a constant at dM. Instead we must appeal to the general solution (2.5), from 
which we infer the asymptotic behavior of a <5£ that preserves the boundary conditions. 
Recall that the leading asymptotic behavior of £ is determined by £o = we®, since w — ► oo 
at X — ► oo. We assume that the boundary conditions are preserved by variations whose 
asymptotic behavior is similar to the sub-leading term in (2.5). This means that we must 
consider the change in the action due to 

5£ = 8Me m (2.19) 

where SM is an arbitrary infinitesimal constant. The corresponding change in the action 
is non-zero and independent of the regulator 

81= J drSM ^ . (2.20) 

This shows that solutions of (2.1) and (2.2) do not extremize the action (1.1) against 
generic variations St; that preserve the boundary conditions on £. In appendix B we repeat 
this analysis in terms of a new metric variable £ that obeys the standard Dirichlet condition 
£ = 1 at dM. Then, in appendix C, we show that the result (2.20) does not depend on 
the diagonal gauge used in (2.3). 

This explains why one cannot drop the regulator-dependent term in (2.15) and assign 
a meaning to the part that remains. The semi-classical approximation of the path integral 
is only well-defined if SI = for the full class of field variations encountered in the path 
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integral. With this in mind, we look for an improved action that is both finite and stationary 
for solutions of the classical equations of motion. 



2.4 The Hamilton- Jacobi Boundary Term 

In the previous section we found that the variational properties of the action are not 
compatible with the semi-classical approximation of the path integral. The same problem 
was encountered in [10], where the authors studied BH solutions of the tree-level equations 
for type OA string theory in two dimensions. In this section we generalize the techniques 
used in [10] to arbitrary theories of dilaton gravity. The result is a new action T that is 
related to / by 

I = T + I CT - (2.21) 
The term Iqt is a boundary integral of the form 

I CT = [ dxL CT {l,X) . (2.22) 
JdM 

The integrand may depend on the fields at the boundary, and their derivatives along the 
boundary, but not their normal derivatives. This ensures that T and I yield the same 
equations of motion. The term Iqt is then defined so that solutions of these equations 
are stationary points of T, with respect to any variation of the fields that preserves the 
boundary conditions. This construction is similar to the boundary counterterm technique 
used in the AdS/CFT correspondence. Although this approach was originally devised to 
cancel divergences in the action [21-24], it was later shown that the boundary counterterms 
follow from the requirement that the action admits a well-defined variational problem with 
asymptotically AdS boundary conditions [19]. A similar modification of the Einstein- 
Hilbert action was proposed in [8,25,26], where the authors considered space-times with 
asymptotically flat boundary conditions. 

Following the derivation in [10], the boundary counterterm is identified with a solution 
of the Hamilton-Jacobi equation for the on-shell action. This equation can be derived by 
starting with the Hamiltonian constraint that follows from the action (1.1) 

H = 2 ir x lab 7r ab + 2 U(X) (j ab vr ab ) 2 + V(X) = . (2.23) 

Recall that the momenta, defined in (2.17), appear as boundary terms in the variation 
of the action. When 51 is evaluated for a solution of the equations of motion the bulk 
term vanishes, leaving only the boundary term. Thus, the momenta can be expressed as 
functional derivatives of the on-shell action 



v^r hab ^ £=oJ x sfj sx 



£=0 



(2.24) 



Rewriting the Hamiltonian constraint (2.23) in terms of these expressions gives a non-linear 
functional differential equation for the on-shell action which is difficult to solve. But in our 
case symmetry arguments reduce the problem to a linear differential equation for a function 
of X, as we now demonstrate. The solution Iqt should be invariant under diffeomorphisms 



- 9 - 



of dA4, so we express its integrand as Lct(i,X) = ^7 Cct(i, X). Because there are no 
other intrinsic, diffeomorphism covariant scalars constructed from 7 in one dimension, it 
follows that Cct can depend only on X and its derivatives along the boundary. But terms 
involving derivatives can be ignored, because the boundary is an isosurface of X. Thus, 
the scalar density Cct depends only on X and the boundary counterterm takes the form 

Ict= / d T ^C CT {X) . (2.25) 

JdM 

The functional derivatives (2.24) applied to this integral give 'momenta' that we denote 

p ab = 1 7 ab £ct{x) px = dxCcT (X) (2.26) 

to distinguish them from n ab and ttx ■ Applying these expressions in (2.23), the Hamiltonian 
constraint becomes a linear differential equation for C CT , 

l - d x (Cct(X) 2 ) + i U{X) Cct(X) 2 + V(X) = . (2.27) 
The general solution of this equation is 

C CT {X) = -y/e-QW {w{X) + c) (2.28) 

where Q{X) and w(X) were defined in (2.6) and (2.7), and c is an arbitrary constant. 
Note that the Hamilton-Jacobi equation determines (Cct) 2 , and the negative root has 
been chosen in (2.28). With the definition (2.21), this root leads to an action T with the 
correct properties. 

In many models there are natural arguments for setting the constant c to zero. For 
instance, in models associated with two-dimensional string theories the action / is invariant 
under Buscher duality [27]. The action V is expected to possess the same invariances as 
/, and this is only the case if c = 0, as pointed out in [10]. More generally, setting c = 
guarantees that Cct is independent of the scaling ambiguity associated with the constant 
Qo that appears in (2.6). Therefore, we set c = and obtain the boundary counterterm 

I C T = - [ dx^Jw(X)e-Q( x ) . (2.29) 
JdM 

This leads to the 'improved' action 

T = I+ [ dx^y\/w(X)e-Q( x ) . (2.30) 
JdM 

2.5 Path Integral Based on the Improved Action 

Based on the analysis of the previous subsections, we conclude that the Euclidean path 
integral for two-dimensional dilaton gravity should weight field configurations using the 
exponential of the improved action, 

Z = f $>g$>X exp(-^T[g,X] \ , (2.31) 



- 10 - 



rather than (1.2). In this section we justify this claim by establishing two important 
properties of T. First we show that solutions of (2.1) and (2.2) are stationary points of 
(2.30) with respect to variations of the fields that preserve the boundary conditions. This 
means that the semi-classical approximation of (2.31) is dominated by solutions of the 
classical equations of motion (2.1) and (2.2). Second, we evaluate the action for the BH 
space-times discussed in section 2.2 and show that the result is finite. 
The variation of the action F is 



5F 



d x yfg 



M 



S^Sguv + SxSX 



+ I dx^ (ir ab - p a0 ) 5 lab + {ti x ~Px)SX 

JdM 



(2.32) 



This differs from (2.16) by boundary terms coming from the variation of (2.29). We can 
repeat the analysis of section 2.3, evaluating 5T on-shell and checking that the boundary 
term vanishes for appropriate variations of the metric and dilaton. The terms involving 
TT ab and ttx were evaluated in (2.18). Using (2.26), the contributions due to Iqt are 



ab 



V 



lj ab e-^VUX), 



Px = \e-^U{X)^UX)- l -^L. (2.33) 



As before, the boundary term is evaluated at the cut-off X = X, 



reg 



6F 



reg 



r/r 



-Q(X) 



+ ( U(X) e"«W (l - J |« ) - W '(X) (l 



(2.34) 



This result simplifies a great deal when the regulator is moved far into the asymptotic 
region of the space-time. The function w(X) becomes very large in this limit, and the 
coefficients of J£ and 6X can be expanded in powers of 1/w. The leading order terms in 
this expansion are 



8T 



reg 



M 
~2 



dr 



1 



^ - U(X) + 



w\x) 

w(X) 



5X 



+ ... 



(2.35) 



The variation of the metric 5£ behaves as (2.19), so the first term clearly vanishes in the 
X rC g — ► cxd limit. However, the analysis of (2.35) is much simpler if we use the change 
of variables = e 2a(x) defined in (B.l), with a = w(X)e Q ^ as in (B.10). Then 
the variation of F is given by (2.34) with all un-hatted quantities replaced by their hatted 
equivalents. Using (B.12), the leading terms in the large-u; expansion of 5F are 



lim 8F 



reg 



M 

2 x v , 



lim 



r/r 



5i + M5 



w 



(2.36) 



The new variable £ obeys the standard Dirichlet condition ^ = 1 at the boundary, which 
means that 5£ vanishes as X reg — > oo. The second term must also vanish in this limit since 



- 11 - 



we assume that w —* oo asymptotically. Thus we have established 

lim 5r reg = 0. (2.37) 

We conclude that solutions of the classical equations of motion are stationary points of T 
with respect to all variations of the fields that preserve the boundary conditions. 

Now we evaluate the action for the BH solution. In section 2.2 we found that the 
regulated on-shell action I reg contains a term that diverges in the X rcg — > oo limit. Including 
the contribution from the boundary counterterm, the improved action is 



r reg = p ( w(X ieg ) J 1 - ™ - w(X ies ) + 2M-2 ttX h T ) . (2.38) 

y y W{Ji re g) J 

The regulator can be safely removed to obtain a finite result 

lim r rcg = [3 (M — 2ttXh T) . (2.39) 

-X"reg — >00 

Notably, this differs by an amount (3 M from the part of (2.15) that was independent of 
the regulator. 



3. Thermodynamics in the Canonical Ensemble 

Motivated by York's path integral formulation of gravitational thermodynamics [11], we 
now introduce an upper bound X c on the dilaton. This constrains the dilaton to a 'cavity' 
X < X c whose wall is the dilaton isosurface X = X c . The cavity is in contact with a thermal 
reservoir that maintains boundary conditions at X c 2 . The semi-classical approximation 
of the corresponding path integral yields the partition function in the canonical ensemble. 
The consistency of this interpretation rests upon the condition of thermodynamic stability, 
which we discuss in subsections 3.5 and 3.10. 



3.1 Dilaton Charge and Free Energy 

In standard thermodynamics the canonical ensemble fixes the temperature and volume of 
a system. The local temperature T c at the cavity wall is fixed by requiring the proper 
periodicity of the Euclidean time to be (3 C . But the volume of the cavity at constant 
Euclidean time is not a well-defined quantity, due to the possibility of a horizon Xh < X c . 
In higher dimensions this issue is addressed by holding the area of the cavity wall fixed, 
instead [11]. We essentially do the same, but in an indirect way. In two-dimensional dilaton 
gravity any sufficiently regular function of the dilaton can be used to construct a conserved 
current. As a result, there are an infinite number of conserved charges associated with the 

2 Note that X c has a clear physical interpretation, unlike the regulator X Icg introduced in the previous 
section. The regulator appeared as an ingredient of the limiting procedure used to evaluate boundary 
quantities. As such, the limit X Icg — > oo was always assumed to be the final step of any calculation. This 
is not the case with X c , which may remain finite. 
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dilaton. Following Gibbons and Perry [12], we choose a particular conserved charge that 
we denote by D{X) and hold this quantity fixed at X c . One convenient choice is 

D C = X C , (3.1) 

which plays a role similar to the area of the cavity wall in higher dimensions. This choice 
can be regarded as "natural" as it is precisely the quantity that couples to the Ricci scalar 
and the extrinsic curvature in the action. Another "natural" possibility for the dilaton 
charge, 

D c = e- Q ^UX c ), (3.2) 

coincides with the solution of the Hamilton-Jacobi equation (2.29). We shall discuss in 
detail how various thermodynamical properties depend on the choice of dilaton charge. 

Fixing D c and (5 C specifies the boundary conditions for the path integral. The classical 
solutions consistent with these boundary conditions correspond to real, non-negative values 
of M that satisfy the condition 

P c = ^Z(D C ,M)P(M) , (3.3) 

where f3(M) is the periodicity of the Euclidean time given in (2.12). In general, there may 
be zero, one, or multiple classical solutions that satisfy (3.3). Issues related to the existence 
of multiple solutions will be discussed in subsection 3.10. For the next several subsections, 
we assume that we are working with a single solution. In that case the Helmholtz free 
energy of the solution is related to the on-shell action T c by 

F C (T C ,D C ) = T C T C (T C ,D C ) . (3.4) 

From (2.14), (2.38) and (3.4) we obtain 

F c = -2tt X h T c + e" Qc (^^0 - ■ (3.5) 

It is worth mentioning that the free energy for the ground state solution vanishes if either 
T = or Xh = 0. This statement does not depend on the location of the cavity wall. 

The interpretation of (3.4) as the Helmholtz free energy assumes thermodynamic sta- 
bility. The consistency of this assumption will be addressed in subsections 3.5 and 3.10. In 
most models thermodynamic stability requires a finite cavity. However, we are often inter- 
ested in the properties of an isolated BH without an enclosing cavity. This is accomplished 
by taking the limit X c — > oo, while simultaneously changing f3 c according to (3.3). This 
limit, if it exists, can be thought of as removing the cavity wall to infinity while adiabati- 
cally relaxing the temperature to preserve M. Motivated by curiosity, we shall sometimes 
apply this limit to the result of a calculation without worrying about whether or not it is 
consistent with thermodynamic stability (usually it is not). We emphasize that, because 
the on-shell action is finite in this limit, any singular behavior in F c as X c — > oo is due to 
the Tolman factor in T c , as evident from (3.4). 

Another issue related to stability concerns the possibility of decay into a field configu- 
ration with a smaller free energy. The BH solution is unstable against decay into any field 



- 13 - 



configuration (that obeys the boundary conditions) with a free energy F g smaller than F c 
from (3.5). We shall discuss this issue in more detail later on. For now, we simply point 
out that there are usually relevant field configurations with a non-positive free energy, so 
a minimal condition for stability of the BH solution is F c < 0. It implies 



2M - 27rX h T < w c ^1 - yjl - . (3.6) 

This condition is always fulfilled for boundary conditions that admit solutions with horizons 
close to the cavity wall: w c = 2M + e. As we shall see in later parts of this section, the 
inequality (3.6) is usually not fulfilled for solutions (3.3) with 2M <C w c . 

Keeping in mind the caveats described above, the collection of formulas in appendix 
D can be applied to the Helmholtz free energy (3.5) to calculate all thermodynamical 
quantities. We start with the most prominent one. 

3.2 Entropy and Area 

The entropy is given by 

(3.7) 



o 9Fq 



The variation of F c at fixed D c can be expressed as 



dF c 



= -27rX h dT c + [ -L - ^ ] dM . (3.8) 



The last term, which was obtained using (D.l) and (D.3), vanishes by the defining relation 
(2.14). Thus, we find the well-known result [10, 12,28-30] 

S = 2irX h . (3.9) 

The entropy is independent of X c , and depends only on the value of the dilaton at the hori- 
zon. As expected entropy is exclusively a property of the horizon. Furthermore, this result 
is universal and applies to any model (1.1) regardless of the potentials U and V. Restoring 
the two-dimensional Newton's constant G2 the result (3.9) becomes S = Xh/(4G2)- It 
is useful to introduce the notion of area at this point in order to compare with the stan- 
dard Bekenstein-Hawking formula. The area of a sphere of radius r embedded in d spatial 
dimensions is given by A d = 2ir d / 2 r d - 1 /T(d/2). In the limit d — > 1 the area becomes inde- 
pendent of the radius and the "sphere" consists of two disjoint points. It is suggestive to 
restrict to a connected component and therefore to associate only one of these two points 
with the horizon. Thus we define = A±/2 = 1. The effective Newton coupling in a 
scalar-tensor theory is given by G e s = G2/X, where X is either approximately constant 
or evaluated at some characteristic scale. The only scale available here is provided by the 
horizon, so we define G e g = GijXh. Then (3.9) is equivalent to the Bekenstein-Hawking 
relation 
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so that the entropy is one quarter of the horizon area in 'effective Planck units'. In section 
6, when we consider the reduction of a spherically symmetric BH from d+1 > 2 dimensions 
down to two dimensions, we shall find that the conserved charge D c always corresponds to 
the higher-dimensional area, in d + 1 dimensional Planck units, of the cavity surrounding 
the BH. This supports our earlier statement about "naturalness" of the choice (3.1). 

3.3 Dilaton Chemical Potential and Surface Pressure 

The chemical potential associated with the dilaton charge (3.1) is given by 

dF c 



dD c 

The variation of F c at fixed T c can be expressed as 



(3.11) 



dF c 



2vT c dX h + I e'^ f^-^)- U c e'^ - y/g) dD c . (3.12) 



Using (D.l), (D.2) and (D.5) we find 

1 „ . a ( ,r 



We have defined the dilaton chemical potential with a minus sign in (3.11) so that it follows 
the same conventions as a pressure in standard thermodynamics. This analogy becomes 
precise in several examples where D c plays the role of the area of a cavity surrounding the 
horizon. 

The dilaton chemical potential (3.13) exhibits some interesting properties: It diverges 
with the Tolman factor as the horizon is approached. For the ground state solution M = 
it vanishes. Finally, it transforms covariantly with the dilaton charge D C (D C ), viz., tp c = 
ip c dD c / dD c . With the choice (3.2) the dilaton chemical potential simplifies to 




= V\It- + (1 -v)\It- - 1 (3-14) 
where 

V ■= : W 1 , ■ (3-15) 

1 - U c w c /w' c 

This definition is particularly useful if rj is constant, such as for models with [7 = 0, 
Minkowskian ground state models and the so-called a6-family. We shall discuss them in 
detail in section 4. 

It is sometimes of interest to determine the values of the dilaton charge for which the 
dilaton chemical potential (3.13) vanishes. This can be thought of as extremizing the free 
energy with respect to the dilaton charge. The dilaton chemical potential vanishes if either 
of the following conditions is satisfied 

■ — / — i — w 1 i — 

tc = Vlo, Vlc = — frV^o- (3.i6) 
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The first possibility can only be realized asymptotically Thus, one natural value of the 
cut-off is the asymptotic region, which might have been anticipated. The second possibility 
arises only if — 1 < w'J (w c U c ) < 0. In the latter case, the system is only stable against 
changes in D c if the free energy is a minimum. This is equivalent to positivity of the 
isothermal compressibility (3.38) discussed below. 

3.4 Internal Energy and First Law 

All calculations so far have referred to the canonical ensemble, with the temperature T c 
and dilaton charge D c at the cavity wall held fixed. We can now perform a Legendre 
transformation to obtain the internal energy E c as a function of the entropy and the 
dilaton charge 

E C (S, D c ) = F C (T C , D c ) +T C S . (3.17) 
From (3.5) and (3.9) we obtain 



The last inequality holds because we assume M, w c > 0. The internal energy E c agrees 
with the proper energy obtained from the quasi-local stress tensor of Brown and York [31]. 
This stress tensor is defined as 

T ab ._ J_ sr_ 

Using (2.32) obtains 

T ab = 2 fab _ p ab\ 



Contracting T ab with two copies of the unit-norm Killing vector u a = 6£ gives the proper 
(Euclidean) energy density. This can be read off from (2.17) and (2.33) 

T ab u a u b = e-^) (V£o " Vt) ■ (3.21) 



This is precisely the internal energy given in (3.18). Notice that it is not the same as the 
conserved charge associated with the Killing vector d T , 



Q 9t = lim y/Z c T ab u a u b = M, (3.22) 

X c — >oo 

where is the lapse function 3 in the 1+1 decomposition of the metric (2.3). The definition 
of the conserved charge Qq t involves the limit X c — ► oo, which should be thought of as the 
two-dimensional analog of the limiting procedure used to integrate over a sphere at spatial 
infinity in higher dimensions. The asymptotic limit of the internal energy is proportional 
to the conserved charge M, rescaled by the Tolman factor 

M 

lim E c = lim — = . (3.23) 



3 Here 'lapse function' refers to the standard temporal evolution, to be distinguished from the lapse 
function of the radial evolution introduced in appendix C. 
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At finite values of X c the result (3.18) can be used to express the conserved quantity M as 
a function of the internal energy 



M= VloE c --^(yToE c ) 2 . (3.24) 
2w c 

This has a simple interpretation in space-times where the Killing norm approaches a 
constant, which we may set to unity without loss of generality. In that case M coincides 
with the ADM mass, which is related to the internal energy in the region X < X c by 

M = E C -^. (3.25) 
2w c 

The ADM mass contains two contributions: the first one is just the internal energy E c and 
the second one is the gravitational binding energy associated with collecting the energy E c 
in a region determined by X c . Consistently, the second term is absent asymptotically. 

Using (D.2) and the results for entropy (3.9) and dilaton chemical potential (3.13) one 
can show that the internal energy obeys the first law of thermodynamics 

dE c = T c dS- i> c dD c . (3.26) 

This quasilocal form of the first law holds regardless of the particular model, the location 
X c of the cavity wall, or choice of the dilaton charge. In this regard it contains a great 
deal more information about the system than the microcanonical expression dM = TdS. 



3.5 Specific Heat at Constant Dilaton Charge and Fluctuations 

Specific heat at constant dilaton charge is given by 



C D 



dE c . 



D c 



dS 



(3.27) 



D c 



Using (3.18), (D.l) and (D.3) we find 



C D 



4tv w' h (w c - 2M) 



2w£(w c - 2M) + w' h 



! 2 



(3.28) 



If w'h then the asymptotic limit of Dp reduces to the microcanonical specific heat 



lim C D = 2ir ^ 



KdMJ 
cPS 
dM- 



(3.29) 



We have implicitly assumed in the previous subsections that the specific heat is positive. 
This assumption will be examined in more detail in subsection 3.10, but we can offer a 
partial justification by considering the near horizon limit of (3.28). If the classical solution 



4 A quasilocal form of the first law was obtained in [32], using background subtraction to address the 
divergences in the action and conserved quantities. 
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determined by (3.3) satisfies w c = Wh + e, with < e < 1, then the specific heat simplifies 
to 

Cd = | • (3.30) 

Remarkably, it is always positive and vanishes as the cut-off approaches the horizon. This 
statement is model independent. 

If the specific heat is positive and the temperature is large enough then one can take 
into account the effect of thermal fluctuations of the internal energy. The temperature is 
'large enough' if quantum fluctuations are considerably smaller than thermal fluctuations. 
Quantum fluctuations are typically of the order A q E c ~ 1, whereas thermal fluctuations, 
A t E c = a/ (E%) — (E c ) 2 = \[C~£)T C ^ depend on specific heat at constant dilaton charge and 
on temperature. Thermal fluctuations dominate over quantum fluctuations if C£>T C 2 S> 1. 
The canonical entropy, which takes into account thermal fluctuations, is given by the well- 
known result 

S fluct = S+UnC D T* + ... (3.31) 
Despite the fact that T c — > oo at the horizon, the fluctuations remain finite there 

lim C D Tl = e~ Qh T. (3.32) 

For some models (including the Schwarzschild BH) (3.32) does not satisfy CoT 2 3> 1, even 
for large BH masses. In these cases thermal fluctuations do not dominate over quantum 
fluctuations near the horizon. 

3.6 Free Enthalpy and Gibbs-Duhem Violation 

With the exception of the dilaton chemical potential, the thermodynamical quantities ob- 
tained in the previous subsections do not depend on the specific choice of the dilaton 
charge. This is not the case for quantities calculated in this and the following subsection. 
Quantities without a bar on top either refer to the dilaton charge (3.1) or are independent 
from it; quantities with a bar refer to the dilaton charge (3.2). 
We start with the free enthalpy (Gibbs function), 

G C (T C , = F(T C , D c ) + ij c D c = E C (S, D c ) - T C S + ^ C D C . (3.33) 

Even though we have Legendre-transformed the internal energy with respect to the quanti- 
ties that usually are considered extensive (entropy and "volume", viz., dilaton charge) the 
result is, in general, not zero. This can be thought of as a 'violation' of the standard form of 
the Gibbs-Duhem relation 5 . There are, however, a few notable exceptions which we shall 
discuss in sections 4 and 5. This shows that extensitivity properties are not necessarily 
as in standard thermodynamics. The physical reason underlying these observations is the 
existence of gravitational binding energy, as in (3.25). Notably, free enthalpy is finite at 
the horizon, despite the fact that free energy diverges there: 



G h = 4?= f 2M - TS + -!- MSU t 



(3.34) 

D c =X h 



For a discussion of alternative forms of this relation relevant for BHs, see [33]. 
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For comparison it is instructive to calculate 



(2r]M-TS) , 



(3.35) 



with r] as denned in (3.15). For constant 77 the quantity \^G C exclusively depends on 
properties of the horizon and the Gibbs-Duhem relation holds if 2rjM = TS. The result 
(3.35) is not only considerably simpler than the expression for G c (which we omitted to 
present) but it is also divergent at the horizon because of the Tolman factor. This is a clear 
demonstration of the fact that free enthalpy crucially depends on the choice of the dilaton 
charge. 

The enthalpy 



H C (S, = E C (S, D c ) + i; c D c = G C (T C , V> c ) + T C S 
can be calculated from free enthalpy (3.35) with the result 



2rj 



M 



(3.36) 



(3.37) 



It is positive if 77 > and also scales with the Tolman factor. The quantity H c is much 
more complicated so we do not present it here. 

3.7 Miscellaneous Thermodynamical Quantities 

The isothermal compressibility 

1 dD c 



kt = — 



D c dif) c 



(3.38) 



requires the calculation of dip c / dD c \j> c performed in appendix D, cf. (D.10). Although 
the dilaton chemical potential diverges at the horizon isothermal compressibility in general 
remains finite there, 



lim kt 



2tt 



(2M) 3 / 2 e^/ 2 



ui e — >2M 



S w"M - 4vr 2 r 2 - AirTMUh - 2M 2 U' h + M 2 U? 



(3.39) 



We do not include the result for kt here, but it may be deduced from the results for Cd 
and C^. 



Specific heat at constant dilaton chemical potential is defined by 



C,/ 



dH, 



dT, 



dS 

Ofr. 



(3.40) 



Using (D.ll) the difference between the specific heats, 

C^-Cd = kt D c T c ( 



(3.41) 
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turns out to be positive if and only if isothermal compressibility is positive. The explicit 
calculation of is somewhat lengthy. A simpler result is obtained starting with the 
dilaton charge (3.2): 



w' 



= 2tt -£ . (3.42) 



w h 



This coincides with the asymptotic limit (3.29) and with the microcanonical specific heat. 
The difference between the specific heats is 



/ 2 



- C D = 2tt -4 • — = ^ . (3.43) 



For generic thermodynamical systems mechanical stability requires positive isothermal 
compressibility or, equivalently [cf. (3.41)], > Co > 0. However, as we have demon- 
strated the quantity depends crucially on the definition of the dilaton charge. Thus, 
unless there is a good physical reason to prefer one particular dilaton charge over all other 
possible choices, we cannot provide an unambiguous answer to the question of mechanical 
stability. It is intriguing that, for the dilaton charge (3.2), thermodynamic stability implies 
mechanical stability: Rx > 0. This is because in a region where both Cp and are 
positive their difference C$ — Co is also positive, as evident from (3.43). 

Additional thermodynamical quantities can be deduced from the ones we have pre- 
sented here. One can basically exploit the results of standard thermodynamics by identify- 
ing the dilaton charge with "volume" , and the dilaton chemical potential with "pressure" . 
A few examples will be provided in section 4. We conclude with a supplementary observa- 
tion regarding an interesting geometric description of the thermodynamics of equilibrium 
systems due to Ruppeiner [34]. The main idea is to introduce a metric on the thermody- 
namic state space, 

ds 2 R := Gf J dX I dX J , G u = -djdjSiX) , X 1 = {E C ,D C ,...), (3.44) 

and to relate its geometric properties to thermodynamic quantities. The entropy S here 
is expressed as a function of extensive variables X 1 , which in our case comprise internal 
energy (3.18), the dilaton charge (3.1), and any additional extensive quantities associated 
with generalizations of the metric-dilaton system. This formalism was applied to various 
BH systems and led to some intriguing results [35]. It is worth mentioning that even in 
the absence of J7(l)-charges the Ruppeiner metric is two-dimensional, which differs from 
the microcanonical analysis in [36], where the Ruppeiner metric was two-dimensional for 
charged BHs only. 

3.8 Charged Black Holes 

Gauge fields in two dimensions have no propagating physical degrees of freedom [37], which 
makes the generalization of our results to charged BHs straightforward. We consider the 
addition of a single abelian gauge field, but the extension to several abelian or non-abelian 
gauge fields is straightforward. The improved action for dilaton gravity coupled to a gauge 
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field is 

r = - \ i d 2 x^g [XR- U(X) (VX) 2 - 2 V(X)] + f, d 2 x^g f{X) F^F^ 
1 JM JM 

- I dx^XK+ [ dxV7e~ Q(X VeoW • (3.45) 
JdM JdM 

The abelian field strength F^ = S/ )M A U — V ' yA^ is coupled to the dilaton field via the func- 
tion f{X). If this function is constant then we say that the gauge field is minimally cou- 
pled, and non-minimally coupled otherwise. We emphasize that the gauge field cannot 
contribute to the boundary counterterm because we require gauge- and diffeomorphism- 
invariance. The only gauge-invariant scalar that is linear in A^ or its derivatives comes 
from the contraction of the field strength with the e-tensor: F fJiU e^ 1 ' . However, this term 
contains a derivative normal to the boundary and is not allowed in a counter-term of the 
form (2.22). All other gauge-invariant scalars constructed from A^, like F^F^, necessarily 
contain such derivatives and are also ruled out. 

The solution of the Maxwell equation V M (/(X) F» v ) = is 

F^ v = - (3 46) 

4/(X) ' { ' 

where q is the electric charge, and the factor of \ has been introduced for convenience. 
When working in the diagonal gauge (2.3) we use the convention e rr = —e TT = +1. The 
equations of motion for the metric and dilaton can be solved in a manner similar to un- 
charged BHs. The metric and dilaton are again given by (2.3) and (2.4), but the square of 
the Killing norm is now 

£(X) = e Qi - x) (w{X) - 2M + ^q 2 h{X)^j (3.47) 

The 'ground state' is defined as the solution with vanishing mass and charge, M = q = 0, 
leading to the same Killing norm i/£o as in (2.8). The function h{X) is defined as 

(X „ Q{X) 

h(X) := J dX , (3.48) 

where the integral is evaluated at X and the integration constant is absorbed into the 
definition of w. The result (3.47) could also have been obtained from our earlier solution 
for the uncharged BH, by replacing V(X) in (2.1) and (2.2) with an effective potential of 
the form 

V cS (X) = V(X) - . (3.49) 

We shall return to this point at the end of this subsection. 
To obtain the on-shell action we use 

ld 2 x^fF^F llu = 2 [d 2 x y /gV„(fF^A„)-2 fd 2 x^g A„V ^fF^) (3.50) 
Jm jm jm 
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and note that the second term vanishes on-shell. The first term, a total derivative, can also 
be interpreted as one-dimensional Chern-Simons terms with support at both the boundary 
and the horizon 



x c 

= -| (CS^Xj-CS^Xh] 

Xh 



2 Jd*x^-gV^fF^A v ) = -| J A T dr 

(3 (A T (X C ) - A T (X h ) ) . (3.51) 



IM 

= _Q 
2 

It is convenient to fix the gauge so that A r = 0. In that case the r component of the gauge 
field is given by 

A T (X) = -| (h(X) - h(X h )) + A T (X h ) . (3.52) 

This is just the potential difference between the cavity wall and the horizon, plus an 
arbitrary constant. The proper electrostatic potential between the cavity wall and the 
horizon is related to A T by 

, , AJX) - A T (X h ) 

* {x) := vm ■ ( ' 

Note that it contains a Tolman factor, like the proper temperature. 

We are now ready to present the generalization of (2.38) to the case of charged BHs, 



T c = p c [e-^ { - VCcJ - 2*X h T c -q<f> c ) . (3.54) 
As pointed out in [10] the logarithm of the partition function, 

Y c (T c ,D C) $ c ) = T c T e = -2nX h T c + e- Q ° (Vto ~ Vtc) ~ q®c , (3.55) 



is not the Helmholtz free energy. Rather, Y c is the Legendre transform of F c with respect 
to the proper electrostatic potential $ c . This is because the Maxwell equation is obtained 
by varying the action with respect to the gauge field, so it is the gauge field (and not the 
field strength) that obeys Dirichlet conditions at dA4. The on-shell action is therefore a 
function of A^ evaluated at the boundary. The derivative of Y c with respect to $ c , holding 
T c and D c fixed, gives the conserved electric charge 



= q . (3.56) 

D C ,T C 



The Legendre transformation 

F C (T C , D c , q) = y c (T c , D c , $ c ) + q $ c (3.57) 

leads to the same expressions for the Helmholtz free energy (3.5) and consequently also for 
entropy (3.9) and internal energy (3.18). 

For minimally coupled Maxwell fields the electrostatic potential $ oc (X c — Xh) is 
proportional to the separation between the cavity wall and the horizon. If the Killing norm 
asymptotes to unity then the proper electrostatic potential $ c diverges in the X c — > oo 
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limit. As a result, the improved action for these models diverges in the limit where the 
cavity wall is removed from the system. We emphasize that this divergence reflects the 
physics of the linear electrostatic potential and therefore cannot be removed from the 
action. The same rule that was mentioned at the end of subsection 2.2 applies here, as 
well. One does not obtain the correct improved action by simply dropping divergent terms. 
In the present case, the action depends on the potential difference between the wall and 
the horizon, which is gauge invariant. Dropping the divergent term from the action leads 
to a violation of gauge invariance. An example of a charged BH with minimal coupling 
is provided by 2D type OA string theory [10], where e~® = w = XX and / = ira! . As 
an example of non-minimal coupling we mention the Reissner-Nordstrom BH, e~® = w = 
2\J~X and f(X) = X (with convenient normalizations for and /). Our formula (3.48) 
gives h = —\/y/X = — 1/r in this case, which is the correct Coulomb-potential. Thus, the 
typical linear (confining) behavior of the electrostatic potential in two dimensions can be 
modified by a non-minimal coupling f{X). 

Maxwell fields open up an interesting possibility. Consider a theory without Maxwell 
fields whose potential V takes the form 

V(X) = V{X) - . (3.58) 

We can also think of this potential as the result of integrating out a gauge field with 
coupling function f(X), in a theory with dilaton potential V{X). We simply re-interpret 
(3.58) as (3.49). These two points of view will generally lead to different conclusions about 
the thermodynamics. In particular, the "natural" choice of ground state is different in the 
two formulations. If we treat the system as dilaton gravity with a Maxwell field then the 
natural ground state is given by M = q = 0. But in the original description, with no 
Maxwell fields, the ground state is associated with M = 0. For instance, the ground state 
M = q = is Minkowski space for the Reissner-Nordstrom BH, whereas in the alternative 
description without Maxwell field the natural ground state M = is a BPS solution. It is 
straightforward to generalize these considerations to potentials of the form 

oo 

V= X nX n/a , (3.59) 

n=— oo 

where a is some natural number. This applies to almost all models in the literature. Now we 
may introduce Maxwell fields to convert all but one of the terms in the generalized Laurent 
series into Maxwell terms. This leads to a rather simple thermodynamical description and 
may be useful in various contexts. We shall provide an example in 6.3 when we consider 
the toroidal reduction of the BTZ BH. 

3.9 Conformal Properties and Microcanonical Thermodynamics 

In appendix B we discuss some consequences of dilaton dependent Weyl rescalings (B.l). 
If interpreted as a conformal transformation, 6 rather than as a change of variables, a new 

6 We use the term 'conformal transformation' in this section, as commonly done in the literature, though 
(B.l) should really be referred to as a Weyl transformation. 
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action (1.1) is obtained with potentials given by (B.3). We say that these models are 
conformally related to each other and sometimes refer to them as representing different 
'conformal frames' within an equivalence class of models. We consider now the transfor- 
mation behavior of thermodynamical quantities. To this end it is very helpful to note 
(B.4) that w is conformally invariant. Therefore, quantities which depend only on w and 
its derivatives are conformally invariant, but quantities that depend on Q are not. 

The conformally invariant thermodynamical quantities comprise the on-shell action T c 
(2.38), surface gravity T (2.13), mass parameter M (2.11) and entropy S (3.9). It is inter- 
esting to note that the specific heat at constant dilaton charge Co (3.28) is conformally 
invariant, despite the fact that both E c and T c depend on the conformal frame. This is 
because the temperature T c (2.14), free energy F c (3.5) and internal energy E c (3.18) are 
conformally covariant, with a conformal weight of —1 due to the Tolman factor 1/y^. 
With this convention the metric is defined to have weight +2. Quantities that are neither 
invariant nor covariant are the dilaton chemical potential ip c (3.14), free enthalpy G c (3.35), 
enthalpy H c (3.37) and isothermal compressibility kt (3.38). These statements are inde- 
pendent from the choice of the dilaton charge. An interesting issue arises for the specific 
heat at constant dilaton chemical potential: while the result for (3.42) is conformally 
invariant, this is not the case for with a generic choice of dilaton charge D c . This 
provides another rationale for working with the charge D c , as defined in (3.2). 

As a simplification one can consider microcanonical thermodynamics (or "horizon ther- 
modynamics") without referring to an actual observer [29]. If one is interested only in the 
microcanonical entropy S, surface gravity T, and their relation to the mass parameter M, 
then one can use a conformal transformation to set U = 0. In this scenario the two relevant 
formulas, besides the result for entropy (3.9), are the Smarr formula for two-dimensional 
dilaton gravity, 

2 

and the microcanonical first law 

dM = TdS. (3.61) 

It is important to realize that the microcanonical thermodynamics is not equivalent to 
canonical thermodynamics in a limit where the dilaton charge approaches some specific 
value, with rare exceptions. The two pictures are equivalent only if there is a value of 
the dilaton charge where the dilaton chemical potential vanishes, ip c = 0, the proper 
temperature coincides the with Hawking temperature, T = T c , and the canonical ensemble 
is well-defined. In that case the canonical first law (3.26) coincides with the microcanonical 
one (3.61), and the mass M coincides with the internal energy E c . The first two conditions 
are satisfied in the asymptotic region of models with a Minkowskian ground state. But the 
third condition is typically violated: in these models the canonical ensemble quite often 
is not defined in the X c — ► oo limit. We show now why this is the case and under which 
conditions the cut-off can be removed. 

3.10 Stability and Phase Transitions 

The most important issue related to stability in the canonical ensemble is the effect of ther- 



mal fluctuations. If the specific heat at constant dilaton charge is negative, Cp < 0, then 
thermal fluctuations in the internal energy have an imaginary component and the system 
is unstable. From the point of view of the path integral, Cd < implies that the Gaus- 
sian integral around a particular extremum of the action does not converge. Fortunately, 
the assumptions w c > Wh and w' h > guarantee that there is an open region around the 
horizon where the specific heat (3.30) is positive. This result is model-independent. Thus, 
if the horizon is located sufficiently close to the cavity wall then the canonical ensemble is 
always well-defined. Of course, a quantitative definition of 'sufficiently close' will depend 
on the features of the particular model. 

In most models a classical solution is only stable when the cavity wall is located at an 
X c that is less than some finite, critical value of the dilaton. Consider a BH associated with 
a solution M of the condition (3.3). We assume that this system has a positive specific heat 
at constant dilaton charge. Now we want to change the boundary conditions X c and (3 C 
while holding M fixed. As described at the beginning of section 3, this is accomplished by 
moving X c away from the horizon while simultaneously relaxing (5 C according to (3.3). The 
equation (3.28) allows us to make a quantitative statement about the point at which this 
process destabilizes the system. Suppose that w is a monotonic function. Then a solution 
X CT it > Xh of 

w' 2 

w(X crit ) = 2M - -\ (3.62) 
2w h 

represents a critical value of the dilaton where the specific heat diverges and changes sign. 
A well-defined canonical ensemble containing the classical solution exists for any value 
Xh < X c < Acrit, but not for X c > X cr it- A similar argument can be made for w that are 
not monotonic. Notice that in either case the specific heat approaches zero at X c = Xh- 
This reflects the fact that the system is only defined for a cavity wall that is external to 
the horizon. 

In general, the condition (3.3) identifies multiple classical solutions consistent with the 
boundary conditions [11]. We refer to each solution as a 'branch' of (3.3), denoted by 
Mj(A c , P c ). The thermodynamic stability of each branch can be determined by looking for 
a solution X l crit > Xh of (3.62). There are three possibilities: 

1. There is no solution A* rit > Xh- This implies that Co is positive, and the branch is 
thermodynamically stable. 

2. There is a solution A* rit > X c > Xh- Again, the specific heat is positive and the 
branch is thermodynamically stable. 

3. There is a solution such that X c > A* rit > Xh- In this case the specific heat is 
negative, and the branch is not thermodynamically stable. 

In the third case the solution Mi{X c , f3 c ) may decay into another field configuration with 
lower free energy via a Hawking-Page phase transition [38]. Because these considerations 
involve only w the emergence of a phase transition happens in any conformal frame. 

We are now able to address the issue of the X c — > oo limit in more detail. Given a 
branch Mj of (3.3) that is thermodynamically stable for some initial values of the boundary 
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conditions, we are interested in the process described by 



X c ^oo lim ^/t(D c ,Mi)/3(Mi) (3.63) 

X c — >oo 

while holding Mj fixed. This removes the cavity wall, leaving an isolated BH. The process 
is only defined if (5 C remains real and non-negative for all values of X c . Furthermore, it may 
be defined on some branches but not others. Finally, even if the process (3.63) is defined 
its endpoint may not be thermodynamically stable. The existence of a thermodynamically 
stable limit (3.63) is far from being a generic feature in the models we consider. For 
instance, the Schwarzschild BH has limx c -»oo Cd < 0. In models where the endpoint 
might be stable, the equation (3.62) can be used to group solutions into two sectors. If 
(3.62) admits a finite solution X CT it > Xh for a particular value of M, then that solution is 
unstable. But if there is no X cr it > Xh for that value of M then the solution is stable. In 
higher dimensions the best known example of this phenomenon is the Schwarzschild- AdS 
solution, where the presence of a cosmological constant stabilizes very massive BHs without 
the need for an external reservoir. An analogous two-dimensional dilaton gravity model 
has wsAdS = 2y/X(l + X/i 2 ) with X = r 2 . For large enough masses, such that Xh > £ 2 /3, 
the quantity w'^ is positive and the system is stable, while for Xh < £ 2 /3 the quantity 
w'h becomes negative and the system is unstable. This can be translated into a critical 
temperature, where the specific heat diverges. It signals a phase transition [38] between a 
BH phase for large masses (small temperatures) and a thermal AdS phase for small masses 
(large temperatures). We shall encounter similar examples in subsection 6.2. 

3.11 Tunneling and Constant Dilaton Vacua 

In addition to the generic family of solutions (2.3), certain models may admit isolated 
solutions known as constant dilaton vacua (CDV). Namely, for each solution X = Xq of 
V(X) = there is a CDV with a line element (2.3) whose Killing norm \/f is given by 

£ = c + ar-^\r 2 . (3.64) 
Here c and a are constants of motion. The constant A is determined by 

r cdv = _q2^ = x = 2V'(X ) . (3.65) 

Because they have a different number of constants of motion CDVs and the generic solutions 
belong to different 'superselection sectors'. Therefore there is no perturbative way in 
which a generic solution could decay into a CDV or vice versa. This is consistent with 
the result [39] that fully supersymmetric solutions necessarily are CDVs, while half BPS 
solutions can never be CDVs. If a supersymmetric CDV exists then the natural ground 
state solution in the generic sector is a half BPS solution, and both kinds of solutions are 
stable. However, once matter is switched on there could be non-perturbative interactions 
which allow tunneling from one sector into the other. We would like to determine whether 
such a tunneling between a generic solution and a CDV is possible, and if so, which of the 
two configurations is favored. 
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We immediately encounter an obstacle: CDVs are isolated in the sense that for an 
open interval of boundary conditions, X c £ (X < ,X > ), there is only a countable number 
of values of the dilaton consistent with the equations of motion. Typically this number is 
either zero or of order unity. Thus, at first glance CDVs are irrelevant because they seem 
to be a set of measure zero in the partition function (1.2). Only for infinite finetuning, 
X c = Xq, there may be a contribution from CDVs. In order to decide whether or not CDVs 
may be disregarded we consider an open interval of boundary conditions for the dilaton and 
the metric. This interval may be arbitrarily small, but it should include exactly one value 
X c = Xq that leads to a CDV. 7 In order to get the number of classical solutions determined 
by this set of boundary data we take the product of the subset of initial conditions consistent 
with the classical equations of motion times the space of continuous constants of motion 
that may be adjusted freely. We call the dimension of this space "weight" and want to 
compare the weight of generic solutions with the weight of CDVs. Let us work in fixed 
physical units. Then, for generic solutions (2.3) there is only a discrete set of values for 
M consistent with the boundary conditions on X c and £ c . So the weight is 2 (from the 
set of boundary conditions consistent with the equations of motion) plus (from the space 
of continuous constants of motion that may be adjusted freely) and thus equals to 2. For 
CDVs the space of allowed boundary conditions is just one-dimensional since the dilaton is 
fixed to a certain value, but now there are two constants of motion, c and a. Therefore, the 
space of constants of motion that may be adjusted freely is one-dimensional. We conclude 
that also for CDVs the weight is 2 and they should not be regarded as a set of measure 
zero. Below we use again sharp boundary conditions and assume X c = Xq for simplicity. 

Now that we have convinced ourselves that CDVs are not a set of measure zero we 
address the question whether tunneling from a generic solution into a CDV is a favorable 
process or not. To this end we compare the generic result for the 'improved' on-shell action 
(2.38) with the on-shell action for a CDV, 

t cdv = _ 2 ^Xqx + ^e-Q^^MXo) . (3.66) 

The first term in (3.66) contains the Euler-characteristic x °f the manifold M. The last 
term is the universal counterterm (2.29) evaluated at X = Xq. It has the same structure 
and, for the same set of boundary data, also the same magnitude as the corresponding 
term in the 'improved' on-shell action for a generic solution. The periodicity in Euclidean 
time is well-defined and positive provided the inequality a 2 + 2cA > holds, which we 
assume henceforth. Minkowski or AdS vacua require a non-vanishing Rindler acceleration 
a while de Sitter solutions (cA > 0) are compatible with a = 0. We conclude that tunneling 
from a generic solution into a CDV is favored if the inequality 

2w(X ) - 2w(X h ) 

Xq X ~ X h > (3.67) 

w h 

holds. Here X^ < Xq is the dilaton evaluated at the horizon of the generic solution. 
Since the right hand side is non-negative the Euler characteristic must be positive for this 

7 In fact, in any experimental situation this smearing of boundary data is physically more reasonable 
than an infinite finetuning. 
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inequality to hold. Therefore we restrict the remaining discussion to the topologies of either 
a sphere (x = 2) or a disk (x = 1). The latter case is consistent with a BH interpretation, 
while the former case may arise for de Sitter. As an example for a model with enhanced 
tunneling to its CDV solution we mention w(X) = c 3 + (X — c) 3 . The CDV solution 
solution for this model is Xq = c > 0, and the inequality (3.67) holds for any value of M 
consistent with Xq > X^. For generic models there are two interesting limits. If Xq is 
close to the horizon, Xq = Xh + e, then for x = 2 the inequality (3.67) holds and tunneling 
to a CDV is enhanced, whereas for x = 1 tunneling is suppressed. If Xq 3> 1 then (3.67) 
simplifies to Xqx > 2w{Xq) /w' h + . . . Since we demand that w — > oo as X — > oo this 
inequality establishes a dichotomy in the space of dilaton gravity theories: if the function 
w diverges asymptotically faster (slower) than linearly tunneling into a CDV is suppressed 
(enhanced), regardless of whether x = 1 or x = 2. 

In summary, a transition between generic solutions and CDVs is possible only through 
tunneling. If the inequality (3.67) holds tunneling from a generic solution into a CDV is 
enhanced, otherwise it is suppressed. We remark that such a tunneling is far from being 
a generic feature: first of all, many models do not exhibit CDVs, and second, many of the 
models that do allow for CDV solutions violate the inequality (3.67) for all possible values 
of Xh < Xq. In fact, none of the models that we are going to discuss explicitly in sections 
4.2, 5 and 6 allows for an enhanced tunneling into a CDV. 



4. Two Dimensional Examples 

4.1 Minkowski Ground State Models 

If the relation 

e Q w = 1 (4.1) 

holds then £ = 1 and the ground state solution (2.8) is (upon Wick rotation back to 
Lorentzian signature) Minkowski space-time. These models are of particular interest as 
they contain asymptotically flat solutions like the Schwarzschild and Witten BHs. As we 
shall see (4.1) leads to considerable simplifications. 
The free energy 

pMGS = E MGS _ ST 



Y 2M 



and internal energy 



Ef GS = w c [l-Jl-— I (4.3) 



2M\ 



Wr 



have a finite limit if the cut-off is removed, lim tUc _» 00 F^ GS = M—TS and lim lt , c _^ tX3 E^ IGS = 
M. The asymptotic version of the first law (3.26) simplifies to the microcanonical dM = 
TdS. The chemical potential conjugate to the dilaton charge D^ GS = w c (3.2), 

TMrq 1 / 2M\ 1/2 1 / 2AA~ 1/2 . . 
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is always non-negative and vanishes for the ground state solution only. It is useful to 
notice that the quantity defined in (3.15) becomes a universal constant, r/ = 1/2. For the 
formulation based upon the dilaton charge (3.1) we get i/j^ gs = 'ifj^ IGS w' c , which is also 
non-negative provided w f c ^ 0. Asymptotically (4.4) implies lim u?c — > oc 

as = M 2 /(2w 2 ). 

From (3.35) we obtain the free enthalpy 

qmgs = M-TS 

2M 



The relevant quantities for isothermal compressibility simplify to f^ GS = w", g^ IGS 
-w"j2 - w' c N c /2, h^ GS = -w"j2 + w' c N c /2 and N c = (2Mw' c w'/ l )/[w c {w' h 2 + 2w£(w c - 
2M))\. Evaluation at the horizon yields lim^—^A/ n^ IGS = 2ir/(Sw'^). In the alternative 
formulation the result is 

^<os _ w .rm . ^-4M+<y< . (4 . 6) 



W r 



AM 2 



The zeros X c > of kt m coincide with the poles of C|f (3.28). We do not investigate 
these models in more detail here but shall encounter them again in this and the next section. 

4.2 Models in the a-b Family 

The so-called aft-family of models [40] , 

U = -~, V = -\x^. (4.7) 

includes the dimensional reduction of spherically symmetric BHs in d + 1 space-time di- 
mensions [a = (d- 2)/{d - l),b = -l/(d - 1)], the Witten BH (a = 1,6 = 0) [1], the 
Jackiw-Teitelboim (JT) model (a = 0,6 = 1) [41,42], and many others. In particular, the 
Schwarzschild BH arises for a = —6 = 1/2. Minkowskian ground state models obey the 
relation a = 1 + 6. Models with a = 1 — 6 have an (A)dS ground state and models with 
6 = a Rindler ground state. Even more general models often are approximated quite well 
by the monomial potentials of the ab family. Therefore our discussion here will reveal some 
generic features of the thermodynamics of dilaton gravity in two dimensions. For the sake 
of definiteness, and to comply with our working assumption (2.10), we assume 6 > —1 and 
B > 0. In addition, we reduce clutter by imposing B = 6 + 1, without loss of generality. 
We fix the constants in (2.6) and (2.7) so that 

w = X b+1 , e Q = X~ a . (4.8) 

With these assumptions the BH horizon is located at 

X h = (2M) 1 /( fc+1 ) . (4.9) 

Note that a negative value of M would result in a naked singularity at X = if a < 2. 
Surface gravity yields 

T = b -±±(2M) b l^ , (4.10) 



-29- 



which leads to a mass-to-temperature relation M oc T l+l l b . The Killing norm squared and 
asymptotic Killing norm squared are given by 

Co = & - 2MX~ a , £o = X b - a+1 > H c . (4.11) 

If we demand / for X c — > oo the inequality 6 — a + 1 > must be met. 

4.2.1 Thermodynamics of the a-b Family 

Entropy (3.9) reads 

5 = 2vr(2M) 1 /( fc+1 ) . (4.12) 

The free energy (3.5) simplifies to 



For b < free energy changes its sign at a particular value of the cut-off, 

1 + 6 
" 1 - 



(4.13) 



£o , (4.14) 



and the asymptotic free energy is positive. This equation always has exactly one solution 
under the assumptions — 1 < b < and M > 0. For Schwarzschild this happens at 
y/X = 9M/4. With the dilaton charge D c = ^•( a + 6 + 1 )/ 2 ^he dilaton chemical potential 
(3.14) simplifies to 

i, c = h \ l M+ ; M-i. (4.i5) 

It vanishes at a particular value of the dilaton given by 

nr l + b nr 



(4.16) 



This equation has a solution if < b + 1 < a. For Minkowskian ground state models 
the dilaton chemical potential vanishes for X c — > oo only. For (A)dS ground state models 
the dilaton chemical potential vanishes precisely at the value of the cut-off that leads to 
vanishing free energy. Internal energy (3.18) is given by 

E c = x( c a+b+l V 2 ^1 - yJl-2MXc 1 ~ b ^j ■ (4.17) 



The specific heat at constant dilaton charge (3.28) simplifies to 

X b+1 - 2M 
bX b+1 - (6 - 1)M 

It changes sign at X CTlt if the inequality 



CD = S _ rb+ c 1 - — . (4.18) 



b _ 1 \W+b) 

X crit = ( M — - — > X h (4.19) 
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holds, which is possible for — 1 < b < only. Thermal fluctuations near the horizon are 
determined by 

lim C D Tf. oc M( a+fe )/( 6+1 ) , (4.20) 

w c ->2M 

with a proportionality constant of order of unity. lib > —a {b < —a) thermal (quantum) 
fluctuations dominate for large masses. If a + b = both kinds of fluctuations are relevant. 
The only Minkowskian ground state model with this property is the Schwarzschild BH. 
Specific heat at constant dilaton chemical potential (3.42), 

CV = f , (4.21) 
is positive for b > and coincides with the X c — > oo limit of (4.18). The free enthalpy 

G c = (l-a-b) b + ] ■ ^= (4.22) 
a + b + 1 y/£ c 

obeys the Gibbs-Duehm relation G c = for (A)dS ground state models only. To get the 
isothermal compressibility (3.38) we calculate (D.6): f c = |(a + 6 + l)(a + 6 — l)Xc a+b 3 ^ 2 , 
g c = |Xr 2 (l - a - X C N C ), and h c = \{b + l)**- 1 ^ + X C N C ), with 

b 2aM - (a - b - l)X} +h /A . 

N c = T-r tjj— ■ (4-23 

2X C {b-l)M-bXl +h V ' 

For any model without (A)dS ground state (a ^ 1 — b) the poles of N c and Cd coincide 
with each other. For (A)dS ground state models N c = b/X c and therefore kt vanishes. 

For Minkowskian ground state models we can discuss the life-time of a BH without 
specifying explicitly the degrees of freedom that carry the Hawking quanta. Let us assume 
that the two-dimensional Stefan-Boltzmann law, asymptotic flux oc T 2 , is valid. Then the 
mass loss per time is determined by dM/dt oc —T 2 . The proportionality constant depends 
on the matter content and on the physical units. With (4.10) it is straightforward to 
calculate the time scale At for the evaporation process 



At oc M^-^+V = M d l^ 



(4.24) 



The second expression in (4.24) refers to the Schwarzschild BH in d + 1 space-time dimen- 
sions. The result (4.24) provides a good estimate of the time-scale associated with BH 
evaporation, as long as the semi-classical approximation remains valid throughout most of 
the process. Typically, this requires that the BH initially have large M. 

4.2.2 Equation of State and Scaling Properties 

For some purposes it is useful to extract an equation of state ijj c (T c , D c ) from the formulas 
above. It can be derived from (4.15) upon expressing the right hand side as a function of 
D c and T c . To this end one has to find M or T as a function of D c and T c using (4.10), 
(4.11) and the defining relation (2.14). As we shall demonstrate, for given T c and D c there 
is in general an ambiguity, i.e., the BH mass M or, equivalently, its Hawking temperature 
T does not follow uniquely from specifying these quantities [cf. the discussion around (3.3)]. 
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A simple exception is provided by Rindler ground state models, where T from (4.10) not 
only is unique but actually is constant, thus establishing 

Rindler : ^{T c , D c ) = — J— 4vrT c D^ a)l{1+a) + ~4t t4^ ^ a " 1)/(a+1) - 1 . (4.25) 

a + 1 a + 1 47ri c 

The dilaton chemical potential becomes minimal at T c = ^JaD^ 1 ^ ( ' a+1 ^/(47r). For the 
Witten BH a = 1 the dilaton chemical potential depends on T c only and attains its mini- 
mum in the asymptotic region. Also models with 6 = 1 lead to a unique T and to a simple 
equation of state, 

MT C , AO = — f 1 + 4^ 2 T c 2 A~ 2a/{a+2) V /2 + — ( 1 + 4vr 2 T c 2 A- 2a/(a+2) r 1/2 " 1 ■ 



a + 2 V / a+2 

(4.26) 

For the JT model a = the dilaton chemical potential becomes independent from the 
dilaton charge. The same happens for all (A)dS ground state models 

{A)dS , te) .i±^ + i^M3_, (4 , 7) 

However, the function k(T c ) in general has more than one branch and is known implicitly 
only. For all other models one has to solve 



T = T C £)(l+b-a)/(l+b+a) J t _ £-2(l+6)/(l+6+a) (^L^j (4.28) 

for T in terms of T c and D c . The solution of (4.28) is not unique in general. For instance, 
the Schwarzschild BH in 5 space-time dimensions, a = —2b = 2/3, leads to two branches 
T± = T c (l ± a/1 - 1447r 2 /(T c 4 L> c ) ) 1/2 / v / 2. If we restrict the discussion to an equation of 
state for the asymptotic region D c — > oo it suffices to solve (4.28) perturbatively to leading 
order, in which case there is only one branch. For Minkowski ground state models the 
result is 

j,2+2/fe 

Minkowski ground state : lim {p c (T c , D c ) oc c - — (4.29) 

D c — »oo J-J c 

while for all other models it reads 

Others: lim ^ C (T C , D c ) oc T 6 1+1 / fe A" (1+fe)(a+6 ~ 1)/(fe(a+fe+1)) • (4.30) 

D c —*oo 

The proportionality constants in these equations are numerical coefficients depending on 
a, b. Asymptotically the equation of state for an ideal gas, ip c oc T c /D c , arises in the limit 
b — ► co, only. As a cautionary remark we recall that the asymptotic limit is not accessible 
for b < because the asymptotic specific heat (4.21) is negative. By contrast, the near 
horizon approximation is always accessible. It corresponds to the limit T c — > oo and allows 
to derive a simple equation of state for the branch where T from (4.28) is finite at the 
horizon, 

Near horizon : ^ C (T C , A) ~ ^ U £)(^-b)/(a+b+i) u 31 ) 

Y y ' a + b + l c y ' 
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Table 1: Extensitivity properties based upon the dilaton charge (3.2) 



Extensitivity properties of BHs are quite unusual; generically, the temperature is not 
intensive and depending on the choice of the dilaton charge, which is extensive by defini- 
tion, either entropy is non-extensive or internal energy or both. For (3.1) the dilaton field 
D c = X c is extensive by definition and thus also entropy is extensive. However, generically 
no other quantity is extensive or intensive. For (3.2) all thermodynamical potentials are 
extensive and the dilaton chemical potential becomes intensive. However, generically en- 
tropy is not extensive. We summarize the extensitivity properties with respect to (3.2) of 
all afr-models in table 1. Extensive (intensive) quantities are denoted by "ext." ("int."), 
while quantities which are neither are denoted by " — " . 

Besides extensitivity properties there are additional interesting scaling properties of 
thermodynamical quantities, which are summarized for generic models in table 2. Here is 
an explanation of all abbreviations: "Conf." refers to the conformal weight (0 means the 
quantity is conformally invariant and "?" means it transforms inhomogeneously) ; "Scale" 
refers to the scale ambiguity — ► ce^, w — > cw inherent to the definitions of these 
functions (0 means again that the quantity is independent from such rescalings); similarly, 
"Shift" refers to the shift ambiguity w — > w+wq (quantities which have no well-defined shift- 
weight are denoted by " — "); "Ext." denotes extensitivity properties (1 means extensive 
and means intensive); "Tol." refers to the scaling with powers of the Tolman factor l/\/£e 
close to the horizon (if it is positive the corresponding quantity diverges at the horizon and 
if the entry is underlined the whole cut-off dependence comes from a single Tolman factor 
with corresponding weight); "Gg 1 " determines the scaling with the gravitational coupling 
constant in front of the action; "Sign" is + if the quantity is positive everywhere outside 
the horizon and "?" if the quantity may change its sign; "Sign^s" is the same as "Sign" 
but for the ground state solution (therefore many quantities vanish); the last entry just 
keeps track of quantities which depend in an essential way on the definition of the dilaton 
charge. All scalings act multiplicatively [s(A) ■ s(B) = s(AB)] except for "Shift" which 
acts additively [s(A) + s(B) = s(A + B)]. The bold entries in this table are definitions 
rather than results; e.g. the quantity is the conformal factor in front of the metric and 
thus by definition has a conformal weight of +2. 

Finally, we collect various restrictions on the parameters from thermodynamical con- 
siderations. We recall that we have assumed B > for positive w, M > for non-negative 
mass and b > — 1 for a well-defined (positive) temperature. If we want space-time to be 
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Table 2: Summary of scaling properties 



regular in the asymptotic region we encounter the inequalities a < 2 and a + b < 1. If 
we demand that the dilaton chemical potential be non-negative everywhere then we need 
< a < 1 + b. Finally, if we wish the specific heat to be asymptotically positive the 
inequality b > must hold. These results are summarized in figure 1. The red (dark gray) 
region is excluded because b < —1 there, thus violating (2.10). The orange (medium gray) 
region exhibits a curvature singularity in the "asymptotic" region X c — > oo. In the green 
(light gray) region the dilaton chemical potential can become negative. In the region below 
the a-axis the specific heat is asymptotically negative. The white region above the a-axis is 
free from thermodynamical pathologies regardless of the location of the cavity. The points 
S (W) [JT] correspond to the Schwarzschild BH (Witten BH) [JT model]. The dashed line 
denotes (A)dS ground state models and the dotted line Minkowskian ground state models. 
The a-axis contains all Rindler ground state models. 

Two particular models of this class deserve further study. The Schwarzschild BH will 
be addressed in more detail in section 6 and the Witten BH will be studied in the next 
section. 

5. Examples from String Theory 

In this section we consider two-dimensional BHs that emerge as solutions, either approxi- 
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Figure 1: Summary of the a6- family 



mate or exact, of string theory. In [1] it was shown that the Euclidean BHs studied in [43] 
admit an exact CFT description in terms of an SL(2,R)/U(1) gauged WZW model. In 
the large k limit, where k is the level of the SL(2) current algebra, the corresponding 
background takes the form 

ds 2 =a'k (dy 2 + tanh 2 y d9 2 ) X = X cosh 2 y . (5.1) 

This solution is commonly referred to as the 'Witten BH'. It is related 8 to the exact 
background obtained in [2] 

ds 2 = dx 2 + v 2 / - dr 2 X = X cosh 2 (6x) Jl - p tanh 2 (6x) (5.2) 

1 — p tanh (bx) v 

where the parameters b and p are 

a b = p := — r = - . (5.3) 

k-2 1 + 2a'b 2 k v ' 

The background (5.2) can be expressed in the form (2.3) by means of coordinate transfor- 
mations described in [44]. 

In the following subsections we wish to consider BHs for the full range k G (2, oo) 
allowed by the CFT. However, in order for the background (5.2) to be a solution of string 



8 The precise relationship is as follows: in the limit k — > oo (or, equivalently. p — > 0) upon identifying 
y = bx and r = V a'k 6 the line-element and dilaton (5.1) emerge from (5.2). 
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theory it must satisfy the condition 

Dim - 26 + 6 a 6 2 = . (5.4) 

Because the target space here is two-dimensional, Dim = 2, requiring the correct central 
charge fixes the level at the critical value k C iit = 9/4. Following [45], we vary k by allow- 
ing for additional matter fields that contribute to the total central charge, modifying the 
condition (5.4) so that k / 9/4. 

5.1 Witten Black Hole 

The Witten BH has the most intriguing thermodynamics of all the models we have exam- 
ined so far [10, 12,28]. We focus on the background (5.1) as a solution of the string theory 
(3- functions at lowest order in a' . The (3- functions can be derived from a target-space action 
of the form (1.1) with 

u = ~ v = ~y x (5 - 5) 

where A = 4/vV. With these potentials and convenient choices for the integrations con- 
stants in (2.6) and (2.7) the functions Q and w are given by 

e Q w = l, w = XX. (5.6) 

Because w is linear in X, the surface gravity and asymptotic temperature T = A/(4-7r) do 
not depend on the mass M. As a result, the condition (3.3) identifies a single solution 
consistent with the boundary conditions (3 C and D c . The relation 2M = XXh then implies 
a simple proportionality between mass and entropy 

M = TS. (5.7) 

An immediate consequence is that the Witten BH satisfies the Gibbs-Duhem relationship: 
GY = Gf = 0. 

The two definitions for the dilaton charge (3.1) and (3.2) are both linear in X c for 
the Witten BH. They are related by D c = XD C . Therefore, the dilaton chemical potential 
associated with D c , 

e^ri-^v /2 + ^i-#r /2 -i, (5.8) 



2 V X c J 2 V X c j 

is related to il>Y from (3.13) by the constant rescaling ij^ = X'ljj^. Both the internal 
energy and the Helmholtz free energy 



Ef = XX c \l-\ll-^\ fY = eY-T™S (5.9) 



have finite asymptotic limits, E w — > M and F w — > 0, respectively. The extensitivity 
properties of the Witten BH are comparable to those of the standard thermodynamics of 
ideal gases, as can be seen from table 1 on page 33. 
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There are various reasons to expect that a semi-classical analysis of the Witten BH 
might encounter problems. For instance, the curvature is of order 1/a! near the horizon, 
indicating that corrections have to be taken into account. It is not surprising that this 
system also exhibits some peculiar thermodynamic properties. The most obvious is the 
behavior of the the specific heat, which is given by 

C™ = 4vr(X c - X h ) . (5.10) 

The specific heat is positive for all X c > X^, but it diverges linearly if we try to remove 
the cavity wall to infinity. In addition, the quantities Cw, and kt are not well-defined. 
However, since we know that the Witten BH is related to the exact solution (5.2), it is 
reasonable to expect that these problems might be absent in a more complete treatment. 
As we shall demonstrate in the next subsection, a' corrections have a significant impact on 
the thermodynamics. 

5.2 String Theory Is Its Own Reservoir: Exact String Black Hole 

The discussion of the canonical ensemble in section 3 required the addition of a cavity 
wall at X c that couples the system to a thermal reservoir. In most cases thermodynamic 
stability is only possible for finite values of X c . We have not questioned the nature of 
this reservoir, or the form of the coupling between the two systems. In previous examples 
we simply assumed that it was possible to couple a given model to some other degrees 
of freedom that functioned as the reservoir. String theory is fundamentally different in 
this regard. It is not possible to couple arbitrary degrees of freedom to the theory, or 
place an abrupt cut-off on the space-time fields. Therefore, if we apply our methods to a 
non-compact solution of string theory any cut-offs introduced during a calculation must be 
removed. 

As shown in the previous subsection, the specific heat of the Witten BH diverges as 
the cut-off is removed to infinity. Remarkably, this and other issues are resolved when a 1 
corrections are taken into account. The metric and dilaton for the ESBH were given in 
(5.2). A target space action of the form (1.1) was constructed for the Exact String BH 
in [44] . In our current notation the relevant functions determining the model are given by 9 

e Q w = l, w(X) = 2b (Vp 2 + 1 + l) , (5.11) 

where b was defined in (5.3), and the canonical dilaton X is related to a new field p by 

X = p + arcsinh p . (5-12) 

It will generally be more convenient to work with the field p, rather than X. The relation 
(5.12) between the two fields is globally invertible, and the limits X — > oo and X — > 

9 Winding/momentum mode duality relates the Exact String BHs to a different set of target space 
geometries, namely the Exact String Naked Singularities. They are also described by (5.11) and (5.12), but 
with sign changes in front of the square root and in front of the arcsinh . 
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coincide with the corresponding limits for p. Therefore, p can be employed not only locally 
but globally. Expressed in terms of p, the square of the Killing norm is given by 



£ = 1 -, • (5.13) 

y^n+i 

This vanishes at 

Ph = ^k(k - 2) . (5.14) 

which establishes the location of the horizon in terms of the level k. 

As in [44], the integration constant in w(X) has been chosen assuming a Minkowski 
ground state, i.e. e®w = 1. Comparing the general form of the Killing norm (2.5) with 
(5.13) shows that the mass parameter is then proportional to the level: M MGS = bk. This 
means that the choice of wq in (5.11) identifies the ground state with k = 0. However, the 
CFT description requires k > 2. This can also be seen from the Hawking temperature [45], 



TESBH =it-l < 5 - 15 > 

which is real and positive only if A; > 2. Therefore, we must shift the constant wq so that 
the ground state corresponds to level k = 2. This is accomplished by replacing w with 
w — 46 in (5.11). Notice that this is the largest shift by a negative constant that is consistent 
with w > 0. The result is 

lESBH(v\ 1 1 



e Q^»{X) = ±_ x W ESBH < 



(X) = 2&(vy + 1-1) • (5.16) 



26 y7TT+i 

Comparing the resulting expression for the Killing norm with (5.13) gives the mass param- 
eter in terms of 6 and k 

M ESBH = _ 2 ) _ ( 5 1? ) 

With (5.16) the dilaton charge (3.2) reads 

D ESBH = 2bp c . (5.18) 

For given j3 c and D c the equation (3.3) has a unique solution with level k > 2. In the 
asymptotic region D c — > oo the temperature is bounded from above and below for k G 
(2,oo) according to (5.15). As we shall show in the next subsection, the thermal partition 
function diverges in the k — > oo limit, which identifies the maximal temperature Tu = 
6/(2-7r) as the Hagedorn temperature. As might have been anticipated on general grounds 
the asymptotic temperature has to be positive and below the Hagedorn temperature for 
the existence of the ESBH solution. 

5.2.1 Thermodynamics of the Exact String Black Hole 

With (5.16) as a starting point we can calculate all thermodynamical quantities of interest 
for the ESBH by applying the general results of section 3. The cut-off is always removed at 
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the end of a calculation, so we drop the subscript 'c' from the following expressions. The 
result for the entropy, 

S ESBH = 2qT f y / k ( k _ 2 ) + arcsinh y/k(k - 2)) , (5.19) 



agrees with equation (4.17) found in [44] (In that reference M = k/2.). For k = 9/4 we 
obtain S ESBH = 27r(3/4 + In 2). In the large k limit the entropy S ESBH = 2irk + O(lnfe) 
asymptotes to the result for the Witten BH with logarithmic a' corrections [46]. The free 
energy, 

f esbh = _ b J^Zl arcs i n h y/k(k - 2) (5.20) 

is manifestly non-positive 10 . This would not have been the case with the original choice of 
wq that lead to (5.11). The result (5.20) coincides with (5.7) in [44]. The internal energy 
is given by the mass, E ESBH = M ESBH = b(k - 2). The specific heat 

C D = C i , = 2irk y/k(k - 2) (5.21) 

is positive and finite. Thus the problem encountered for the Witten BH is absent: removing 
the cut-off is possible for the Exact String BH and leads to a positive and finite specific 
heat. It is worth mentioning that the specific heat and entropy both vanish in the limit 
T — > 0, so the Exact String BH is one of the few examples for which the third law of 
thermodynamics holds. 

Even though the Minkowski ground state condition is not satisfied by (5.16), the Killing 
norm still approaches unity as the cut-off is removed. Therefore the asymptotic limit of 
(3.15) is constant, j] ESBH = 1/2. As a consequence the dilaton-chemical potential vanishes 
as X~ 2 in the limit X c — > oo. The Gibbs-Duhem relation 

qESBH _ pESBH _ j^ESBH _ rpESBH gESBH _^ q ^ 22) 

is violated even in the large k limit, 

lim Q ESBH = lim F ESBH = -b\u2k + O ( \ In k ) . (5.23) 

k-^oo k—>oo yk J 

Thus, comparing the thermodynamic behavior of the exact String BH with that of the 
Witten BH demonstrates the importance of a' corrections, even in the large k limit. In 
the limit k — > 2 the Gibbs-Duhem relation, 

lim G ESBH = lim F ESBH = -b{k - 2) + 0{k - 2) 2 , (5.24) 

is restored. 

It is interesting to observe how the free energy scales with the level k in various 
limits. For large values of k the free energy scales logarithmically, while for k close to the 



It is worth mentioning that there is also a CDV solution with X = and AdS geometry with curvature 
proportional to b 2 . It has vanishing free energy. Therefore the Exact String BH is stable against tunneling 
(cf. subsection 3.11). 
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ground state value it is linear in k — 2. This implies corresponding behavior of the thermal 
partition function (2.31). Suppose that we fix the factor in front of the action as a/(2ir), 
where a is some positive parameter, so that the curvature term in the action 1.1 is given 
by — aXR/ {An). Then the thermal partition function in the large k limit is given by 



gESBH 



= (2fcr + ...~ (S KSBH ) +... (5.25) 

fc>2 



In this limit the partition function grows as a positive power of A;. It diverges in the k — > oo 
limit, as the temperature approaches the asymptotic value Th = b/(2-rr). In this sense, Tjj 
represents the Hagedorn temperature. Similarly, for k — > 2 we obtain 



rESBH 



k^2 



Thus, in the weak coupling limit the partition function grows only with a power of entropy, 
while in the strong coupling limit it grows exponentially with entropy. This is consistent 
with fermionization in the strong coupling regime. 11 For arbitrary values of k we obtain 

■gESBH n\ _ e <r arcsinh-y/fc(fc-2) ^ 

Because exp [arcsinh (3/4)] = 2 we obtain the interesting result 

Z ESBH {9/4) = 2 a . (5.28) 

As a caveat, we remind the reader that these results apply to the semi-classical approxima- 
tion of the thermal partition function, based on the action for the ESBH given in [44] . They 
do not take into account contributions from the tachyon that appears in the bosonic string 
spectrum, or from any matter that must be coupled to the ESBH to maintain the correct 
central charge (5.4). Furthermore, while we have focused on the effects of a' corrections, 
we have not taken into account contributions from world-sheet topologies other than the 
leading term. 

5.2.2 Thermodynamical Derivation of the Exact String Black Hole 

It is interesting to note that the potentials (5.11) can be derived purely from thermody- 
namical considerations. Of course, this requires some ad-hoc assumption: we look for a 
model which may be shifted to a Minkowski ground state model and whose BH solutions 
have a specific heat whose asymptotic behavior is 

C D = oc M 2 T (5.29) 



11 In the WZW formulation the stress tensor is given by 

T(z) = -*—r, ab J a J b + -(d^f + bdc. 
k — I 4 

Here rjab is the metric on the algebra sZ(2,M) and the J a are given by equation (2.11) in [2]. For k — > 2 
the first term dominates and the system behaves like a free Fermi gas, cf. e.g. [47]. This suggests that an 
appropriate interpretation of the Exact String BH in that limit is not in terms of some effective geometry 
but rather in terms of free fermions. 
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for all values of M and T. The rationale behind this assumption comes from the con- 
sideration of the weak {k — > oo) and strong (k — ► 2) coupling limit. In the former case 
we have seen that the ESBH is approximated by the Witten BH which has a divergent 
asymptotic specific heat. Taking into account corrections from fluctuations at 1-loop level 
leads to a specific heat Cd oc M 2 T [48]. This is a non-stringy calculation, but as long as we 
are not interested in the precise proportionality factor it should not matter which kind of 
fluctuations we consider perturbatively: thermal fluctuations, a' corrections and quantum 
fluctuations of some additional massless matter fields are expected to yield the same qual- 
itative behavior. In the strong coupling limit the system can be fermionized (cf. footnote 
11). A free Fermi gas has a specific heat of the form C oc T. The proportionality constant 
can be rescaled such that it includes a factor M 2 . Then the weak and strong coupling 
limits qualitatively both yield the same result for the specific heat, which motivates our 
Ansatz (5.29). 

With Cd oc w' h /w'^, T oc w' H and 2M = Wh the most general model (1.1) obeying 
(5.29) has to fulfill the autonomous differential equation 

w" oc • (5.30) 

w 

Without loss of generality we fix the proportionality constant in (5.30) to unity. Its general 
solution is 

dw 

c ±X. (5.31) 



\Jc\- 2/w 

The integration constant c\ may be adjusted by rescaling w, Co and X, so we fix it to unity. 
The integration constant co may be adjusted by shifting the origin of X, so we fix it to 
zero. The sign ambiguity is resolved by requiring X to be positive. Thus we obtain 

y/w(w-2) + In w - 1 + ^w{w - 2) =X. (5.32) 

Defining p = w(w — 2) reveals that (5.32) actually is equivalent to (5.12). Solving w in 
terms of p and taking the branch of positive w leads to 



W (X) = V V + 1 + 1 , X = p + arcsinh p . (5.33) 

This is the same as w(X) in (5.11) for b = 1/2. If we had chosen a different proportionality 
constant in (5.30) we would obtain a different value for b but otherwise the same result for 
w. The function Q is fixed by the Minkowski ground state property (4.1). Thus, the Exact 
String BH in the Minkowski ground state representation (5.11) is the most general two- 
dimensional dilaton gravity model with a Minkowskian ground state that has an asymptotic 
specific heat Cd = oc M 2 T. 

6. Spherically Symmetric Black Holes in Higher Dimensions 

Up to this point we have focused on theories in two dimensions, but the results of section 
3 can also be applied to solutions of gravitational theories in higher dimensions, as long 
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as the on-shell action reduces to the form (2.38) or (3.54). With a few caveats, this leads 
to a simple description of the thermodynamics of a class of spherically symmetric BHs in 
d + 1 > 2 dimensions, with or without a cosmological constant. 

Pure gravity with a cosmological constant in d + 1 > 2 dimensions is described by the 
action 

= [ d d+1 x^i{R d+1 -2k)-—^ ( d d x^T d {K d + ...) . (6.1) 

Subscripts l d+ 1' and '<f have been used to distinguish quantities from their two- and one- 
dimensional analogs. The '. . .' in the boundary integral indicates boundary counterterms, 
whose precise definition depends on the space-times asymptotics [8, 19]. We are interested 
in solutions of this theory that can be expressed as a direct product of the two-dimensional 
metric (2.3) and the metric on a round d — 1 sphere 

ds 2 = i{r) dr 2 + ^dr 2 + (G d+1 )^ ip{r) 2 dfi£_i . (6.2) 

Explicit factors of the Newton's constant have been included in this expression so that the 
field f(r) is dimensionless. Ignoring for a moment the contributions from the boundary 
counterterms, the on-shell reduction of (6.1) on the d— 1 sphere leads to a two-dimensional 
expression of the form (1.1). The dilaton X in this action is related to the scalar field (p by 

X(r) = T G d+1 tp(r) d - 1 T := (6.3) 

8 7T U d+ l 

where A d _\ is the solid angle subtended by the d — 1 sphere, and the constant T has 
been defined for convenience. Thus, the dilaton X(r) is the proper area of a sphere with 
coordinate radius r, in d + 1 dimensional Planck units. The kinetic and potential functions 
for the dilaton are 

V{X) = ~\{d- l)(d - 2) T^i Jffcf + e X . (6.5) 

The last term in V comes from the cosmological constant, which we parameterize in terms 
of a length scale I by 

did -I) , . 

A:=6 e = ±1 '°- ( 6 - 6 ) 

The coefficient of the reduced action explicitly sets the two-dimensional Newton's constant 
to 87rCr2 = 1. From (6.4) and (6.5) we obtain 

i 

e Q(X) = i X3 ^ ; g 7 

d — 1 

= (<i - 1) Jl - 1 T^n J , (6.8) 

where we have fixed the free constants in the definitions (2.6) and (2.7) conveniently. 
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The results of section 3 can be applied to these solutions of the d + 1 dimensional 
theory. This is because the on-shell action, after integrating over the sphere, takes the 
form (2.38). Of course, the contributions from the boundary counterterms, alluded to by 
the '. . .' in (6.1), must be properly accounted for. Otherwise the reduced on-shell action 
would diverge, as in (2.15). This raises an interesting technical problem. The formulation 
of the boundary counterterms in the d + 1 dimensional theory depends on the asymptotics 
of the space-time [8,49]. This would seem to imply that the reduction of theories with 
A = and A / must be studied separately. We take a different approach. Rather 
than working out the reduction of the appropriate boundary counterterms in (6.1), we 
assume that the two-dimensional boundary counterterm derived in subsection 2.4 gives a 
sensible renormalization of the reduced on-shell action. For the solution described above, 
the boundary counterterm is 



I CT = -p w (X) ] jl-^. (6.9) 

This counterterm gives the 'correct' on-shell action in the sense that the analysis of section 
3 recovers the standard results associated with the thermodynamics of BHs whose metrics 
take the form (6.2). An immediate consequence, due to (6.3), is 

S = 2<KX h = . (6.10) 

The entropy of these BHs is always given by one-quarter of the proper area of the horizon, 
in d + 1 dimensional Planck units. 

There is no reason to believe that the counterterm (6.9) can be lifted to higher di- 
mension, to recover the counterterms of the original theory. The combination of the two- 
dimensional metric and scalar field are simply not sufficient to reproduce the appropriate 
functional dependence on the higher- dimensional metric. But despite this caveat, it seems 
remarkable that the quasi-local thermodynamics of spherically symmetric BHs with differ- 
ent space-time asymptotics can be recovered from the reduced on-shell action, renormalized 
by the two-dimensional counterterm (2.29). 

6.1 Schwarzschild Black Hole in d + 1 > 4 Dimensions 

As a first example, consider the on-shell reduction of (6.1) with A = 0. Setting (p ~ r, up 
to factors of the Newton's constant, gives the standard form of the Schwarzschild solution 

2M , 

X(r) = Tr d - 1 (6.12) 
where the mass parameter M is related to the constant M by 

M:=^. (6.13) 
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In section 3 the dilaton charge D c = X c and proper periodicity (3 C were held fixed at r = r c . 
This is precisely the canonical ensemble studied in [31,50], where the BH is enclosed inside 
a cavity of radius r c whose area and proper temperature are held fixed by means of an 
external thermal reservoir. The relation (6.3) identifies the dilaton charge D c = X c as the 
area of the cavity, so the dilaton chemical potential (3.13) is the corresponding pressure 



The internal energy of the system can be expressed as 



(6.14) 



E c = (d-l)Tr c d 



2M 



(6.15) 



The asymptotic limit of this quantity gives the ADM mass M. Furthermore, the internal 
energy satisfies the quasi-local form of the first law (3.26), 



dE c (S,A c ) =T c dS- Pc dA c 



(6.16) 



where A c is the area of the cavity in Planck units (8irGd+i = 1)- Finally, the expression 
(3.28) for the heat capacity recovers the well-known result for the thermodynamic stability 
of the Schwarzschild BH enclosed by an isothermal cavity in d+1 dimensions 



C D = — 4tt M rh — ^ 



2M 



rf 2 -dM 



This quantity is positive for <r c < r max , as indicated in figure 2. 

c D 



2M 



(6.17) 



Figure 2: Heat capacity for the Schwarzschild BH in dimension d + 1 > 4. 



6.2 Asymptotically AdS Space-Times 

Now we turn our attention to the case of non-zero cosmological constant. For definiteness, 
we take the cosmological constant to be negative (e = —1) and consider asymptotically 
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AdS space-times 12 . Setting tp ~ r leads to the Schwarzschild-AdS solution 

r 2 2 /If 

&) = 1 + p~^=2 X = Tr d -K (6.18) 

As in the previous case, holding D c = X c and (3 C fixed at r c corresponds to an AdS 
BH enclosed by an isothermal cavity. The thermodynamics of this system was studied 
in [51,52]. 

Rather than reproduce the full catalog of results regarding the Schwarzschild-AdS 
solution, we focus on a few interesting properties which are qualitatively different than the 
previous example. First, the internal energy (3.18) is given by 




Unlike the result (6.15), the asymptotic limit of the internal energy vanishes due to red- 
shifting 

lim E c = . (6.20) 

r c ^oo 

However, as in section 3.4, the conserved charge is [31] 

Q 9t := lim v^) E c = M . (6.21) 

r c — >oo 

This is the same result obtained in higher dimensions using any one of a number of methods 
[53]. But, notably, it does not agree with the standard AdS/CFT result when the space- 
time dimension d+ 1 is odd. In that case, global AdS (the ground state) has Qq t = M cas , 
where M cas is a Casimir energy associated with the dual field theory on S 1 x S* ' -1 [22,23]. 
Reproducing the Casimir term is possible in the present context: an appropriate shift in wo 
can be used to shift the value of (6.21) by an amount M cas , for all solutions. But from the 
two-dimensional point of view the shift is arbitrary. The counterterm (6.9) is completely 
ignorant of the physics attached to this issue in the higher-dimensional theory. 

Another interesting result is the structure of the heat capacity for an asymptotically 
AdS BH. The result is, in general, very complicated for finite r c [52]. But the presence 
of a cosmological constant can stabilize the BH even when the cavity wall is removed to 
infinity. Using (3.29), the r c — ► oo limit of the heat capacity is 

lim C D = ( d - 1) 2,X h { ^+('-2)^ , (6 22) 

Unlike the asymptotically flat solutions studied in the previous subsection, the heat capacity 
may be positive for a large enough BH. The relation between the heat capacity and the 



12 Asymptotically de Sitter space-times with A > can also be addressed in this framework, but many 
of the working assumptions made in sections 2.1 and 3 must be changed. We shall discuss this briefly in 
section 7. 
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size of the BH is shown in figure 3. The condition Cd > is satisfied if the horizon radius 
is larger than a critical size set by the dimension and the AdS length scale 



d-2 
d 



This corresponds to a critical periodicity for the Euclidean time given by 

y/d{d - 2) 



crit 



2tt£ 



(6.23) 



(6.24) 



The possibility of a phase transition, as studied by Hawking and Page [38], is encountered 
when rh drops below the critical value given by (6.23). This is consistent with our general 
discussion in subsection 3.10. 



d = 3 d 



d=2 





r h =£ 



Figure 3: Heat capacity for the Schwarzschild-AdS BH in d + 1 dimensions, with the cavity wall 
sent to infinity. The horizontal axis is the horizon radius 77,, and the vertical axis is the X c — > 00 
limit of Cd- 



6.3 Reduction of BTZ 

Gravity with a negative cosmological constant in 3 dimensions admits a solution known as 
the BTZ BH [15]. If A = -l/£ 2 then (Euclidean) line element of the BTZ BH is 



ds 2 BTZ = i{r)dr 2 + -^-dr 2 + r 2 (# - ^-dt) 2 (6.25) 

where 



r 2 1 6(7 2 T 2 
£(r) = -.-8G 3 M+ 1—. (6.26) 



If the inequality | J\ < Ml holds then two horizons exist, 
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The normalization chosen in the rest of this work would require G3 = 1/4, but we set 
G3 = 1/8 to obtain the same normalization of the action used in [15,54]. Some useful 
formulas are 



Achucarro and Ortiz performed a Kaluza-Klein reduction 

ds 2 BTZ = g af3 dx a dx f3 + X 2 (d<p + A a dx a f (6.29) 

to two dimensions. The dilaton field is the surface radius, X = r, and the two-dimensional 
line element g a pdx a dx^ corresponds to the rr-part of (6.25). They integrated out the 
gauge field A a and obtained an action of the type (1.1) with [54] 

U A ° = 0, V AO = -§ + ^~ 3 - (6-30) 

We point out an important deficiency of (6.30): the angular momentum J enters here as a 
parameter in the action rather than emerging as a constant of motion. Consequently, the 
thermodynamics is not reproduced correctly from the two-dimensional model (1.1) with 
the potentials (6.30). But there is a simple solution to this problem: we have demonstrated 
in subsection 3.8 how to treat charged BHs and noted that upon introducing Maxwell fields 
a parameter in the action may be converted into a constant of motion, namely a conserved 
charge. In fact, the Kaluza-Klein Ansatz (6.29) already contains such a Maxwell field and 
all one has to do is to refrain from integrating it out. Therefore, we study an action (3.45) 
[multiplied by 2 to get the normalizations above] with the functions 

U BTZ {X) = 0, (6.31) 

V BTZ {X) = ~§, (6-32) 

f BTZ {X) = - l -X\ (6.33) 

The potentials (6.31) and (6.32) are consistent with (6.4) and (6.5), respectively, for d+1 = 
3. The Maxwell field with coupling function (6.33) arises because we are treating a spinning 
BH rather than a spherically symmetric one. The conserved U{\) charge therefore has to 
be identified with the angular momentum, q oc J. The proportionality factor depends on 
the conventions. With our choices the correct relation is q = J/2. For further calculations 
we need the functions 

Q BTZ (X) = , w BTZ (X) = * , h BTZ {X) = A (6.34) 

and note that the ground state solution M = J = is pure AdS3, £0 = X 2 /I 2 . 

We show now that this effective two-dimensional description leads to the correct ther- 
modynamics for the BTZ BH. To this end we apply the general results of subsection 3.8 
to the specific choice (6.34) recalling that q = J/2. The square of the Killing norm 

£(X) = *! - M + ^ (6.35) 
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indeed is equivalent to (6.26) for G3 = 1/8 and X = r. This together with the gauge 
potential 

A r = 0, A T = -^, ( 6 - 36 ) 

plugged into (6.29) allows to recover the BTZ line-element (6.25). With the proper angular 
momentum [cf. (3.53)] 

tt ° " *- m th ~ *?) (6 37) 

the thermodynamic potential (3.55) is given by 

Y(T c ,X c ,n c ) = -T c S + E c -n c J. (6.38) 

Entropy is determined by the value of the dilaton at the outer horizon, Xh = J"+. Taking 
into account that the normalization of Newton's constant differs in this subsection from 
the one used in the rest of this work by a factor of 2 we get from (3.9) 

S = 4vrr+ . (6.39) 

The result (6.39) coincides with the entropy of the BTZ calculated by counting of mi- 
crostates, cf. e.g. [55]. The inverse Euclidean periodicity determines temperature, 



y.2 y-»2 



2M / J 2 



1 " 775*5 • ( 6 - 4 °) 



2nr + £ 2 S V M 2 l 2 

This coincides with the result for the temperature of the BTZ BH [15]. Obviously, in 
the limit \J\/{Mt) — > 1 temperature vanishes, which is consistent with the fact that for 
I J I = Ml the BTZ BH becomes extremal. If we define the quantity O as 

0:= lim JIsIc = 4y ( 6 - 41 ) 
X c -»oo Ar+ 

then we obtain the correct first law relating the conserved charges, 

dM = TdS + ttdJ . (6.42) 

The integrated version 

2M = TS + 2QJ (6.43) 

shows that extensitivity properties differ slightly from standard thermodynamics. However, 
they are as expected for a BH in AdS: if the dilaton (and thus entropy) is extensive then 
the square-root of the mass M is extensive, consistently with the extensitivity properties 
in table 1. This explains the factor of 2 on the left hand side of (6.43). The factor 2 on 
the right hand side is the standard result for spinning BHs. With the results above, the 
expression for internal energy [note again the factor of 2 as compared to (3.18) because of 
the different normalization of the action] , 

/— 2X ( I MP T 2 f 2 \ 
E c (T e ,X e ,J) = 2yfe-2yfe=^[l-Jl-^ + ^\ , (6.44) 
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and surface pressure 



E, 



'c 



J 2 



(6.45) 



•c — 



lead to the quasilocal form of the first law 



dE c = T c dS + Q c dJ — tf) c dX c 



(6.46) 



All our results are compatible with [51] (cf. also [56]). This shows that the BTZ BH thermo- 
dynamics is compatible with a Kaluza-Klein reduction to two dimensions. A generalization 
to a charged BTZ BH is straightforward and requires the addition of another Maxwell field 
with a linear coupling to the dilaton. A different generalization involves the coupling to 
a gravitational Chern-Simons term [57], which yields upon dimensional reduction [58] a 
model of type (3.45) [59] whose entropy recently has been derived in [60] using Wald's 
Noether charge technique [61]. 

7. Discussion 

In this paper we have studied black hole thermodynamics for arbitrary models of dilaton 
gravity in two dimensions. The analysis of these models was complicated by the fact 
that the path integral defined with respect to the action (1.1) does not have a sensible 
semi-classical limit. We resolved this problem by constructing an improved action, 



generalizing the techniques first used in [10] to obtain the boundary counterterm (2.29). 
The counterterm is unique up to a shift ambiguity 13 which we fixed by demanding that it 
be independent under changes of scale associated with constant shifts in Q{X). Thus, one 
should consider (7.1), instead of (1.1), as 'the generic dilaton gravity action'. Following 
the approach of York [11], we proceeded to define the canonical partition function for a 
black hole enclosed in the 'cavity' X < X c . The Helmholtz free energy (3.5) for this system 
is related to the improved action by F c = [3~ l T C: where the proper inverse temperature 
/3" 1 is fixed at the cavity wall X c by an external thermal reservoir. All of the results 
on thermodynamics in section 3 follow from this identification. After deriving various 
thermo dynamical quantities, we generalized our results to include two dimensional Maxwell 
fields and their associated conserved charges. We then examined the conditions under which 
the system is thermodynamically and mechanically stable, and described processes which 
might destabilize the system. Our results were then applied to a wide range of models in 
sections 4-6. This includes backgrounds that appear in two-dimensional solutions of string 
theory (section 5). Finally, we discovered that our results can be used to describe the 
thermodynamics of certain higher dimensional black holes (section 6). 

13 In addition there is a sign ambiguity which is trivially fixed by demanding finiteness of the improved 
on-shell action. 
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An important aspect of our analysis is universality. For instance, we confirmed that 
the entropy (3.9) for a dilaton gravity black hole in two dimensions is essentially given by 
Xfu the dilaton evaluated at the horizon. This result is independent of the details that 
define a particular model. If we define the 'area' of the horizon via the d — ► 1 limit of a d 
dimensional sphere and restore Newton's constant, then this result can be expressed as 



The effective Newton's constant is defined as G c s = GijX^. Since G2 is dimensionless 
there is no notion of Planck length in two dimensions. However, we can think of G e g as 
defining an effective resolution at the horizon, in the sense that decreasing G e s increases the 
number of bits that can be stored by the black hole. Then (7.2) can be interpreted in much 
the same way as the usual relationship between entropy and area in higher dimensions. In 
section 3.4 we obtained a simple expression for the (proper) internal energy of the system 



This agrees with the calculation of the energy using the quasilocal stress tensor defined by 
Brown and York [31]. We were able to show that every dilaton gravity black hole satisfies 
a quasilocal form of the first law of black hole thermodynamics 



This expression holds at arbitrary values of X c , incorporates the dilaton charge D c and its 
chemical potential ijj c , allows for Maxwell fields with proper electrostatic potentials $ c and 
conserved charges q, and accounts for the (non-linear) effects of gravitational binding energy 
present in E c . The canonical first law (7.4) contains considerably more information than 
the microcanonical dM = T dS, which becomes a triviality if expressed as dwh = widX^. 

Despite our best attempts to provide a comprehensive study, there are several direc- 
tions that we have not had a chance to explore. A number of 'working assumptions' were 
made throughout the paper and it is important to understand which ones can be relaxed. 
An example is our decision to focus on models where w — ► 00 as X — > 00. This is a rea- 
sonable assumption, as it seems to be the most common behavior among dilaton gravity 
models. But there are certainly models for which it does not apply. For instance, one 
might consider models where w approaches a constant, or w — > —00, as X — > 00. Our 
approach should still be applicable, though many of the considerations in section 3 must 
be modified. In the first case, where w approaches a constant, we can make a constant 
shift so that w — > as X — > 00. Then one can exploit the duality described in [62], which 
preserves the classical solution for the metric but replaces w with w = 1/w. This implies 
w — > 00 as X — > 00, which can be analyzed under our initial assumptions. The second case, 
w — > — 00 as X — > 00, is perhaps more interesting as it includes the dimensional reduction 
of the static patch of de Sitter space. We now have to consider the region X < XdSi 
where X<is corresponds to the observer-dependent cosmological horizon in the static patch. 
Although we have not performed this analysis, we expect that our approach should repro- 
duce the standard results for the thermodynamics of cosmological horizons. An additional 




(7.3) 



dE c = 



T c dS - ip c dD c + <& c dq . 



(7.4) 
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complication emerges when one considers a second horizon, Xh < X^s, associated with a 
black hole. The analysis of that system is beyond the scope of this paper. 

It would be interesting to study the impact of the boundary counterterm in (7.1) 
on the exact path integral quantization procedure in the first order formulation of dilaton 
gravity, cf. e.g. [3,63]. While the bulk term and the Gibbons-Hawking- York boundary term 
look quite different from (1.1) in the first order formulation, we expect that the boundary 
counterterm (2.29) can be translated easily to its first order counterpart, because it depends 
solely on the dilaton field and the induced volume form at the boundary. 

Another idea that deserves to be studied in more detail was mentioned at the end of 
subsection 3.8. Consider a dilaton gravity model with a potential of the form (3.59). The 
same model can be obtained by integrating out a number of Maxwell fields whose coupling 
functions f{X) are related to terms in the potential. This means that various coupling 
constants in the potential V{X) can be converted to conserved U (1) charges associated with 
Maxwell fields. This technique allowed us to study the thermodynamics of the BTZ BH in 
subsection 6.3. The first law for the rotating, three dimensional solution was recovered by 
'integrating in' a gauge field that Achucarro and Ortiz had integrated out in their toroidal 
reduction [54]. It would be interesting to study the extent to which the thermodynamics 
of higher dimensional solutions can be reproduced in this manner. 

The most pressing generalization of our work involves coupling two dimensional dila- 
ton gravity to matter, such as a scalar field. This is important for several reasons. As 
mentioned in the introduction, we have assumed that these models will be coupled to 
some form of propagating matter, and scalar fields are a prime candidate. We have also 
considered backgrounds that arise in two dimensional solutions of string theory, and any 
further analysis of these models should try to incorporate the tachyon. In general, coupling 
additional fields to the dilaton gravity model will introduce complications at the level of 
the Hamilton-Jacobi equation. However, even if this equation can no longer be solved 
exactly, it can still be addressed perturbatively. This approach has been used successfully 
for gravity coupled to a scalar field in higher dimensions [64-67] . We hope to address this 
in a future publication. 
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A. Conventions, Dimensions, and a Dilaton Gravity Bestiary 



Let Ai be a two-dimensional manifold with boundary dAi. We consider a metric g^ on 
Ai that induces a metric j ab on dAi. The (one-dimensional) boundary indices are useful 
for keeping track of factors of 7 and its inverse. 

The Riemann tensor on (A4,g) and its contractions are defined so that spheres have 
positive curvature and hyperbolic spaces have negative curvature 

R\ P u = dpT^ + T^T x Kp -(u~p) (A.l) 
Rfiu = R X fjiXu ■ (A- 2) 

The extrinsic curvature associated with the embedding of (dM.,~j) in (At,g) is given by 

K ab = l - n^d^ab (A.3) 

where is an outward-pointing unit vector normal to dAi. The trace of the extrinsic 
curvature is K = j ab K ab . Because the boundary is one- dimensional, it follows that 

K ab - lab K = 0. (A.4) 

This implies that Tr ab would vanish if not for the coupling between X and the Ricci scalar 
in (1.1). 

Throughout the paper we work with Euclidean metrics of the form 

ds 2 = £(r) (It 2 + — dr 2 . (A.5) 
£(r) 

Calculating the action and its variation requires the Ricci scalar and the extrinsic curvature 
of the boundary. The Ricci scalar is given by 

d 2 £ 

R = -—^ = - e ~Q[ w " + Uw l + U'(w-2M)] . (A.6) 

The last expression, which follows from (2.4) and (2.5), is coordinate independent and 
relates the Ricci scalar to the dilaton field. For models obeying the Minkowskian ground 
state condition (4.1) it simplifies to R = 2M(w' 2 /w — w"). Boundary quantities are ob- 
tained from the X — > 00 limit of a constant X surface, as described in subsection 2.2. The 
induced metric on the surface is j ab = £,(X), and the trace of its extrinsic curvature is 

^■^..^(-y, (A . 7) 

2 dr 2Vw-2M y ' 

For models obeying the Minkowskian ground state condition (4.1) it simplifies to K = 
MvJ /(wyf£). The Euler characteristic for the space-time is related to (A.6) and (A. 7) by 

4vrx= d 2 x^/gR + 2 / dx^K. (A.8) 
JM JdM 

One can verify that the BH solutions of discussed in 2.2 have the topology of a disk, \ = 1- 
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The dimensions of various quantities are assigned conveniently, but at times differently, 
in the examples we consider. Here we clarify the role of physical dimensions, how to fix 
them, and how to rescale them. Unless stated otherwise we always work in natural units 
h = c = k-Q = 1 and choose the dilaton field X to be dimensionless. This means that the 
potential U{X) is dimensionless, as well. If the coordinates are taken to have the natural 
dimension of length then the relation (2.4) implies that must also have dimensions 
of length. This can be accomplished by defining = Aexp / dXU(X), where A is a 
dimensionful constant with units of length that includes the factor of e^°. In that case, 
both w and M must have units of inverse-length, so that (2.5) is dimensionless. Similarly, 
coordinates with units of length imply that the function V has dimensions of one over 
length squared. Units of length can then be changed by appropriate rescalings of the 
dimensionful constants in V, as well as A and M. Alternatively, if the coordinates are 
taken to be dimensionless then U, e®, V, w, and M are all dimensionless, as well. 

After these general observations we focus on the choices of relevance for our work. With 
one exception (mentioned below) we demand that the dilaton field X be dimensionless. 
In some of our examples we do not exhibit the dimensionful constants, for simplicity. 
Instead we treat all quantities as dimensionless, including the line-element, and restore 
appropriate physical dimensions when necessary. This applies to the discussion of the 
aft-family in subsection (4.2), and to the discussion of the BTZ BH in subsection (6.3) 
if X,£,J,M are considered to be dimensionless. In other examples we keep track of the 
physical dimensions for easier comparison with the literature. This applies to sections 5, 
6.1, and 6.2. Alternatively, in subsection 6.3 the quantities X,£,J can be thought of as 
having dimensions of length, with M dimensionless. 

The analysis of sections 2 and 3 is carried out for a general model of two-dimensional 
dilaton gravity. As a convenience, the properties of several models are summarized in tables 
3 and 4. Table 3 lists models that contain only a metric and a dilaton, while table 4 lists 
models that also contain Maxwell fields. The last three models in table 3 can be converted 
into models belonging to table 4 as explained at the end of subsection 3.8. In most entries 
in these tables we have chosen the scaling and shift ambiguities inherent to Q and w as 
to achieve coincidence with the choices in the main text. In some brief examples in our 
paper the four-dimensional Newton constant is fixed as G4 = 1/2, thus simplifying the first 
entries in both tables, as well as the Schwarzschild-(A)dS entry in table 3. In subsection 
4.2 we set B = b + 1, which simplifies the corresponding entry 'ab- Family' in table 3. 

B. Weyl Rescaled Metric 

In subsections 2.3 and 2.5 we employ a change of variables that involves a Weyl rescaling 
of the metric. Then, in subsection 3.9, we consider the effect of a conformal transformation 
on a model. In this appendix we work out some of the relevant details. 
Consider a new metric related to g^ v by 

9^ = e 2a Wg^ (B.l) 
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Model 


U(X) 


V(X) 


e Q\X) 


w(X) 


Reference 


Schwarzschild 


1 

2X 


1 

2G 4 


/ Ct4 

V 2X 


p2X 
V Gi 


[68] 


Jackiw-Teitelboim 





-AX 


1 


AX 2 


[41] 


Witten BH 


1 

X 


-¥* 


1 
AX 


XX 


[1] 


CGHS 





A 
2 


1 


XX 


[69] 


(A)dS2 Ground State 


a 
X 


-\BX 


1 

X a 


B V"2-a 
2— a 


[70] 


Rindler Ground State 


a 
X 


-\BX a 


1 

X a 


BX 


[71] 


BH Attractor 





B 
2X 


1 


BhxX 


[72] 


a6-Family 


a 
X 


B Y a +b 
2 A 


1 

x a 


B vb+l 
5+1 A 


[40] 


Liouville Gravity 


a 


be aX 


e aX 


2b Ja+a)X 
a+a c 


[73] 


Exact String BH 


[44] 


[44] 


(5.16) 


(5.16), (5.12) 


[2] 


Schwarzschild- (A) dS 


l 

2X 


1 3 y 

2G 4 e 2 ^ 


/ Ga 
V 2X 


(l + ^A) 


[38] 


Katanaev-Volovich 


a 


PX 2 - X 


e aX 


I e aX {X-(5X 2 + 


[74] 














KK Reduced CS 





|X(c-X 2 ) 


1 


\(c-X 2 f-\c 2 


[58, 59] 



Table 3: Summary of Dilaton Gravity models, adapted from 



Model 


U(X) 


V(X) 


f(X) 


e Q(X) 


w(X) 


h(X) 


Reference 


Reissner-Nordstrdm 


l 

2X 


l 

2G 4 


X 


[gI 

V 2X 


Fix 

V Gi 




[75] 


Achucarro-Ortiz 





X 


X 3 
8 


1 


X 2 


4 

X 2 " 


[54] 


2D Type OA/OB 


1 

X 




na' 


1 

AX 


XX 


lnX 


[76] 



Table 4: Dilaton Gravity models with Maxwell fields. 

where cr(X) is an arbitrary function of the dilaton. If we express the action (1.1) in terms 
of this new metric we find 
1 



i[x,g^} 



d xy g 



M 



XR - U(X)(VXy - 2 V(X) 



dxy/^XK. (B.2) 
8M 



This has the same functional form as (1.1), with dilaton potentials U and V that are related 
to U and V by 



U(X) = U{X) - 2a'(X) V(X) = e 2a{x) V(X) 



(B.3) 
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The functions Q{X) and w{X) transform according to 



Q(X) = Q(X) - 2 a(X) w(X) = w(X) . 



(B.4) 



A change of variables with 2a{X) = Q{X) leads to vanishing U. This transformation is 
common in the dilaton gravity literature because it simplifies many calculations. Notice 
that w(X) is invariant under (B.l). Referring to (2.15), which depends only on w(X), X, 
and the parameter M, we confirm that the on-shell action is invariant. This is as it should 
be: the value of the action cannot depend on the variables we use to evaluate it. 

Now we want to examine the variation of the action using the new metric variable. In 
terms of g^ u and a a general variation of g^ u takes the form 

Sg^ = e 2a Sg^ + 2 e 2a g^a'5X . (B.5) 

The change in the action (2.16) due to a small variation of the fields can be written as 



SI[X, g^v 



d x\/g 



M 



S^Sgav + SxSX 



+ 



dxy 7 
8M 



(B.6) 



Both the equations of motion and £x, and the momenta 7r and ttx have the same 
functional form as in the main text, with all un-hatted quantities being replaced by their 
hatted counterparts. The solution of the equations of motion is 



X = X(f) 



£(r) dr 2 + 



1 



dr 2 



d f X 



-Q(X) 



£(f) 

i(X) = e^ x ) (w(X) - 2 M) 



(B.7) 
(B.8) 



This is just the solution from section 2.1, expressed in terms of the new metric and written 
in a coordinate system (r, f) with df = e 2a ^ x ^d r . Evaluating (B.6) for this solution gives 



SI[X,g^} 



£=0 



dx 
8M 



1 



8£ + e 



(B.9) 



We can now choose the function &{X) in (B.l) to simplify the analysis of this expression. 
If we define 

e 2 ^= W (I)eW, (B.10) 



then the (square of the) Killing norm for the new metric is given by 



w(X) 



(B.ll) 



This makes the meaning of Dirichlet boundary conditions unambiguous: £ takes the con- 
stant value £ = 1 at dAi, and any variation <5£ that preserves this condition must vanish 
at dM. The functions Q and U are determined as 



a Q(X) 



1 

MX) 



u 



w'(X) 

~<X) 



(B.12) 
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so that (B.9) yields 



5L, 



dr 



w 8£ 



1 



w J 



(B.13) 



The subscript indicates that this quantity is evaluated using a regulator, as in section 2.3. 
Following the arguments in the main text we ignore the 5X term and concentrate on a 
variation of the metric 5£ = 5 Me® = SM/w. The corresponding change in the action is 



lim SL 



reg 



drbM / 



(B.14) 



which of course agrees with (2.20). Changing variables to a metric that obeys standard 
Dirichlet conditions at dA4 simplifies the analysis of the improved action, as well, and is 
useful in establishing the result (2.37). 

In subsection 3.9 we treat (B.l) as a Weyl transformation, as opposed to a change of 
variable. Then the action (B.2) is interpreted as a functional of X and g^v- This gives a 
new model of type (1.1) with potentials given by (B.3). 



C. Proof of Gauge-Independence 

In this appendix we show that the main results of section 2 are independent of the choice 
of gauge (2.3). We begin by establishing two important properties of solutions of (2.1) and 
(2.2). Consider the vector k^ 1 defined by 

k» = e^ u e Q{x W u X (C.l) 

where e^ x ^> was given in (2.6). The Lie derivative of the dilaton along k^ vanishes by 
construction 

C k X = e Q( - x) e^ u d^XdvX = . (C.2) 

The same result also applies to the metric. It follows from the equation of motion (2.1) 
that the vector k^ satisfies Killing's equation 

C-hQixv = V^ky + Vyk^ = . (C.3) 

Therefore, solutions of (2.1) and (2.2) always possess at least one Killing vector whose 
orbits are isosurfaces of the dilaton. 

The properties of k^ allow us to construct a coordinate system (r, r) where the Killing 
vector is k^d^ = d T , and the metric and dilaton take the form 

ds 2 = N(r) 2 dr 2 + £(r) (dr + N T (r) dr) 2 X = X(r) . (C.4) 

The functions N and N T that appear in the metric are the lapse function and shift 'vector', 
respectively, associated with evolution in the direction orthogonal to k^. Normally, the 
lapse and shift both appear in the action as Lagrange multipliers. In this case, the fact 
that all fields are constant along the direction singled out by the Killing vector implies that 
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the action is independent of N T (r). Furthermore, the coefficient of 5N T in the variation of 
the action vanishes identically. This means that we can freely set N T = in (C.4) with no 
loss of generality. The metric then simplifies to 

ds 2 = N{r) 2 dr 2 + £(r) dr 2 . (C.5) 

One can substitute this form of the metric into the action and recover all of the results of a 
fully covariant analysis. From the point of view of the action, the functions £(r) and iV(r) 
must be treated as independent fields. However, we are always free to present a solution 
in a coordinate system where N(r) 2 = ^(r)" 1 , as in (2.3). 

Now we want to check whether or not the conclusions of subsection 2.3 are sensitive 
to the choice of gauge. Integrating the equations of motion for an arbitrary lapse function, 
we find 

8 r X =N(r) v / £(X) e~ Q{x) 

(.-«). 

The only difference between these expressions and the solution (2.4) and (2.5) is an overall 
factor that appears in d r X. Notably, we obtain the same expression for £(X), concurrent 
with its coordinate independent interpretation as square of the Killing norm. The change 
in the action due to a variation in the metric, with the bulk terms set to zero, is now given 
by 

5I ICR = -- [ dr 9rX 5£(r) . (C.7) 

rCg 2 J iY(r)v^) { ' 

The subscript is a reminder that the integrand is evaluated using a regulator. This result 
can be obtained either by evaluating the covariant expression (2.16) for the solution (C.6), 
or by expressing the action (1.1) in terms of £(r) and N(r) and then varying the fields. 
It is clear from (C.6) that the explicit factor of N(r) in (C.7) cancels against the factor 
in d r X, and we recover the result (2.20). Therefore, the conclusions of subsection 2.3 are 
independent of the choice of gauge. This exercise can be repeated for the variation of the 
action T from (2.30). Regardless of the choice of gauge, ST = for all variations of the 
fields that preserve the boundary conditions. 

D. Formulas 

The variations of X^ and T can be expressed entirely in terms of dM, using the definitions 
of the horizon (2.11) and T (2.13) 

1 1 c)T 

dX h = —— dM dT = —— — — dM . (D.l) 

2ttT 2ttT dX h v ' 

Similarly, the definitions of £ (2.5), Q (2.6), and w (2.7) give 

di c = -2 e Qc dM + (U c Cc + e Qc w' c ) dD c . (D.2) 
When D c is held fixed the variation of £ c reduces to 

dC c =-2e Qc dM. (D.3) 
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A useful consequence is 



die 
dT 



-4ire 



w 



D c w h 

Holding the temperature T c fixed implies a relationship between dM and dD c , which can 
be obtained using the definition (2.14), and the relations (D.l) and (D.2): 



(D.4) 



T c 2 

Using (3.5) and the formulas above we obtain 
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The relations (D.4) and 
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allow to determine 
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dT c 

Another useful formula is 
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